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Abstract. We consider series of the form ^2 a n {n- x}, where n S 1 d and {x} 
is the sawtooth function. They are the natural multivariate extension of Dav- 
enport series. Their global (Sobolev) and pointwise regularity are studied and 
their multifractal properties are derived. Finally, we list some open problems 
which concern the study of these series. 



1. Introduction 

Let [ • J denote integer part and let { • } be the centered sawtooth function denned 

by 

(x-[x\ -i if x$Z 
] else. 

The purpose of this paper is to investigate regularity properties of the multivariate 
functions which are defined by 

\Jx G R d f(x) = ' x }> ( 2 ) 

where n ■ x denotes the standard inner product between the vectors n and x, and 
(a„).„ eZ d is a real valued sequence indexed by the set Z d = Z d \ {0}. With a 
slight abuse, the vectors n for which a n is nonvanishing will be referred to as the 
frequencies of the series. 

In the one- variable case, examples of such functions can be traced back to the Ha- 
bilitationsschrift of Riemann, see [101113]; they were later considered by Hecke |28j . 
and also Hardy, who studied the series 

°° S X 

%(*) = £^- (3) 

n— 1 

It seems however that the general one-dimensional case was first considered only 
in 1937 by H. Davenport in [19l [20]. The first of these papers starts with the 
following remarkable identity, which establishes in all generality the connexion with 
Fourier series: 

> a n {nx\ = y Cm sin(27rmx) with c m = } na n . (4) 

^— ' ^— ! -Km 

n=l m=l »a 

One of the fascinating aspects of these expansions is that their study lies at the 
crossroad between several areas of mathematics. They appear naturally in several 
problems related with analytic number theory; this actually was the motivation of 
H. Davenport for studying them, see also the recent studies by R. de la Breteche 
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and G. Tenenbaum (such as |17| for instance). They were later considered in con- 
nexion with harmonic analysis, see e.g. |35j and references therein where a function 
space point of view is developed, and it is shown in which sense an arbitrary one- 
periodic odd function can be expanded on this system. Convergence properties of 
these series at particular points are related with the Diophantine approximation 
properties of these points, see [17[ [35] . Recently, J. Bremont studied the L 2 and 
almost-sure convergence of these series, see [T5] . The multifractal analysis of these 
functions shows connexions between their pointwise regularity properties and geo- 
metric measure theory, see |35j , and also |47j for an extension of Davenport series 
with translated phases. Note also that examples of Davenport series valued in 
K 2 were proposed by H. Lebesgue as space filling functions; this study was devel- 
oped in |38[ 139] , where the connexions between Davenport series and space filling 
functions are examined. 

In this paper, we shall investigate the multivariate case, which has not been con- 
sidered up to now. Our main motivation is that multivariate Davenport series are 
natural examples of multifractal fields. The recent increase of interest in such fields 
is motivated by the relevance of multifractal analysis techniques in image classi- 
fication, see [U [2]. Indeed, the validation of 2D multifractal analysis algorithms 
requires the introduction and the mathematical study of collections of multifractal 
fields of various kinds. However, very few multivariate multifractal models have 
been studied up to now (see however the PhD thesis of H. Oppenheim [41] for an 
early analysis of a multifractal function of several variables where Diophantine ap- 
proximation properties are involved, and [4] for fields generated by random wavelet 
series). Another case of random fields which have recently been studied are Levy 
fields, which are a natural extension of Levy processes to the multivariate setting; 
their multifractal analysis has recently been performed by the authors, see [24 . The 
scarcity of existing results is partly due to the fact that the derivation of the mul- 
tifractal properties of multivariate functions lies on variants of ubiquity methods 
which can prove much more involved in the multidimensional setting. Therefore, 
extending the collection of available multivariate models, and elucidating their mul- 
tifractal properties is an important issue. An additional motivation of this paper 
is to draw a comparison between the multifractal behavior of Davenport series and 
Levy fields. Indeed, they are both constructed as superpositions of piecewise linear 
functions which display jumps along hyperplanes, the main difference being that 
the locations of these hyperplanes are random in the case of Levy fields, whereas 
they are determined by arithmetic conditions in the case of Davenport series. We 
shall see that the multifractal properties of Davenport series bear similarities with 
those of Levy fields, so that they can be seen as a kind of deterministic counterpart 
of these fields. 

The paper is organized as follows. In Section [2j we establish the relationships 
between the Davenport and Fourier coefficients of Davenport series, in the normally 
convergent case. We shall see that this relationship extends to more general, and 
actually distributional, settings in Section |8l The main purpose of this paper is 
the study of pointwise regularity properties of Davenport series. The key step 
consists in analyzing the locations and magnitudes of the jumps of Davenport series. 
Preliminary results concerning this study are collected in Sections [3] and 2] In 
Section^ an upper bound of the Holder exponent is derived and, as a consequence, 
cases where this exponent vanishes everywhere are worked out. A difficult question 
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(which is far from being closed, even in the one-variable case) is to understand 
when this upper bound is sharp; the purpose of Section |6] is to show that this 
is the case when the frequencies of the Davenport series are sufficiently sparse. 
Implications for multifractal analysis are stated in Section [7] In Section [8] we 
shall consider convergence properties of Davenport series in the Sobolev spaces H s 
for set, especially when the sequence of coefficients does not belong to I 1 ; this 
study will also be the occasion to draw bridges with arithmetic functions in several 
variables, a topic which has been barely scratched until now (see however [TBI |2"5] 
and references therein). Concluding remarks and open problems are collected in 
Section [9] Finally, the proofs of the main results are completed in Sections [10] 
and 1111 This is the occasion for us to exhibit deep connexions with the theory of 
sets of large intersection and with the Duffin-Schaeffer and Catlin conjectures in 
the metric theory of Diophantine approximation, see Section 111.21 

2. Relationships between Davenport and Fourier series 

We start by establishing some conventions. First, note that functions such as ([2]) 
are necessarily odd and Z d -periodic. Since {— x} = — {x}, it follows that the system 
supplied by the {n ■ x}, for n € Z rf , is redundant. The choice made for one- 
dimensional Davenport series is to use only these functions for n > 1, as e.g. in Q 
above. We shall make a different choice in dimension d > 2, which will preserve the 
symmetry of the decomposition. Specifically, we shall keep both functions {n ■ x} 
and {— n ■ x}, and, without loss of generality, we shall assume that the sequence 
(cbn)nez d °f Davenport coefficients is an odd sequence indexed by Z d , which implies 
uniqueness of the decomposition. 

The function spaces that we shall consider are composed of Z d -periodic odd 
functions, and the sequence spaces that we shall consider are composed of odd 
sequences. Therefore, we shall use the following conventions concerning spaces: 
With a slight abuse of notations, £ p will denote the space of odd sequences which 
belong to £ p (Z d ), L 2 is the space of odd locally square- integrable functions which 
are Z d -periodic, and, more generally, if E is a space of functions defined on M. d , we 
shall also denote by E the space of odd functions that belong locally to E and are 
Z d -periodic. 

Another convention concerns divisibility in several dimensions. Let n, m € ZJjf. 
If m = In for some ! e Z*, we say that I and n are divisors of to. The fact that 
the term "divisor" applies without distinction to elements of 7L d and of Z* will not 
create confusions because the context will always be clear. If I = ±1 are the only 
integer divisors of to, we say that to is irreducible; this means that its components 
are coprime. Throughout the paper, N denotes the set of positive integers. 

Finally, the support of a sequence a = {a n ) n ez d is 



Let us now investigate the relationship between Davenport and Fourier scries. To 
this end, let us assume that the series ([2]) converges normally, i.e. that the sequence 
( a n)nGZ d belongs to £ l (convergence properties in different functional settings will 
be investigated in Section [8]). Then, / belongs to L°°, hence to L 2 and the Fourier 
series expansion of / converges in L 2 . Since / is odd, it may be written in the form 



supp (a) = {n e Z rf | a„ ^ 0}. 
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Here, we adopt the same convention for Fourier series as for Davenport series, that 
is, we assume that the expansion is taken on all frequencies of Z d , but that the 
sequence (c m ) meZ d is odd. Since ^ is normally convergent, 




for all m e Zf. A straightforward computation shows that the above integral is 
equal to zero except if n is a divisor of m, in which case there exists an integer 
I £ Z* such that In — m, and the integral is equal to — l/(2irl). As a consequence, 
the Fourier coefficients of / are given by 

(i,Ti)ez* xz d 

ln = m 

Note that without making any assumption on the summability of the sequence 
(a n ) ngZ d, the above formula still enables us to define a sequence {c m )mez d - This 
detour via Fourier series will allow us to study the convergence of the series ^ 
even when (a n ) n Gi d does not belong to i 1 , see Section[S] Indeed, we shall see that, 
in many functional settings, when the associated Fourier series converges, then the 
partial sums of the series ^ n a n {n ■ x} converge to the same limit. 



1 



3. Discontinuities of Davenport series 

Let us consider a bounded function g : R d — > M. By definition, the magnitude of 
the jump of g at any fixed point xq € K d is 

A- g (x ) = limsup <7(x) — liminf g(x] 

X^XO X—f-XQ 

(which may possibly vanish, in which case g is continuous at xq). The magnitude 
of the jumps can also be expressed by means of local oscillations. To be specific, 
let us recall that the oscillation of the function g on a bounded subset f2 of K d is 
defined by 

Osc fl (f2) = sup g(x) — inf g{x). 
leu zen 

Letting B(x, r) denote the open ball with center x and radius r, it is easy to see 
that the magnitude of the jump of the function g at the point xq satisfies 

A s (a; ) = lim Osc g (B(x , r j). 

We shall now determine the set of points at which the Davenport series / defined 
by ([2]) has a discontinuity, and we shall study the magnitude of the corresponding 
jump. We assume in what follows that the sequence (&n) n ez d belongs to t l . 

Given a vector with integer coordinates g € and an integer p 6 Z, let H p q 
denote the hyperplane 

H Ptg = {x€R d \p = q-x}. (6) 

It is clear that multiplying p and the components of q by a common integer value 
leaves the hyperplane unchanged. In order to ensure the uniqueness of the repre- 
sentation, it is sufficient to assume that p and the components of q are coprime and 
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that q belongs to the subset Zl of III formed by the vectors whose first nonvanish- 
ing coordinate is positive. In fact, one easily checks that any hyperplane H Piq may 
be indexed in a unique manner by a pair (p, q) that belongs to 

H d = {(p,q) eZxZ d + \ gcd(p, <?) = !}, 

where gcd(p, q) is the greatest common divisor of the integer p and the components 
of the vector q. Furthermore, let { • denote the restriction of the sawtooth func- 
tion { • } to the open interval (—1/2, 1/2). Then, { • }* is continuous everywhere 
except at the origin: {x}+ makes a jump of size —1 when x crosses zero in the 
upward direction. In addition, {x} is the sum of {x — p}+ over all the integers 
p G Z. Along with the fact that the sequence (a n ) n ez d an d the function { • }* are 
both odd, this enables us to rewrite the definition ([2]) of the Davenport series / in 
the form 

oo 

f( x )= X! fp,i( x ) with fp,q( x ) = 2 ^2 a h{K<l ■ x -P)}*> 
(p,q)en d i=i 

where all the series converge normally. This decomposition enlightens the fact that 
/ is the superposition of a family of functions that are continuous everywhere except 
maybe on a specific hyperplane of the above kind. To be precise, each function 
is continuous everywhere except maybe on the hyperplane H pq , and the fact that 
(p, q) belongs to Hd implies that these hyperplanes are distinct. Moreover, when 
a point x crosses H pqi the real points l(q ■ x — p), for I > 1, all cross zero, so that 
f P ,q{x) makes a jump of magnitude \A q \, where 

oo 

A g = 2j2ai q . (7) 

i=i 

Note that, in the case where the latter sum vanishes, the function is actually 
continuous on the whole space, including the hyperplane H pq . 

The analysis of the discontinuities of the Davenport series / begins with a first 
remark: As the series ((2|) is normally convergent, its sum / is a function in L°°, so 
that the potential discontinuities must have finite magnitude. We shall now show 
that the set of points at which / is not continuous is exactly 

U H v-i- (8) 

First, note that, if a point xq does not belong to the latter set, then the above 
decomposition entails that / is a sum of uniformly convergent series of functions 
that are continuous at Xq, thereby being continuous at xq as well. Conversely, if 
a point Xq belongs to a hyperplane H p q indexed by a pair (p, q) € Hd for which 
A q does not vanish, and to no other hyperplane of that form (which is the case of 
Lebesgue-almost every point of H p ^ q ), then / has a discontinuity at xq of magnitude 
\A q \, that is, 

A f (x ) = \A g \ > 0. 

More generally, suppose that xq belongs to a (possibly infinite) collection of hy- 
perplanes H puqi indexed by pairs {pi,qi) € Hd for which A qi do not vanish. The 
previous case shows that, for any specific value of i, the function / has a discontinu- 
ity of magnitude exactly \A qi \ on a dense set of points of H Pi , qi . Therefore, one can 
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pick a point yi arbitrarily close to xq at which / has a discontinuity of magnitude 
exactly | A qi \ . It follows that 

A/(iro) > max|^ 9i | > 0, 

i 

so that / exhibits a discontinuity at xq- 

Given that \A q \ is the magnitude of the jump of the Davenport series / at 
Lebesgue-almost every point of the hyperplane H pq , we shall call with a slight 
abuse \A q \ the magnitude of the jump of / on H p q . 

We see here a sharp contrast with Fourier series: The series ^ will usually 
exhibit discontinuities no matter how fast the coefficients a n decay. The following 
proposition shows that even more is true: The zero function is the only continuous 
Davenport series. 

Proposition 1. Let f be a Davenport series with coefficients given by a sequence 
a = ( a n) n ez d * n ^ ■ If f i- s a continuous function, then 

Vn e Z d a n = 0. 

Proof. The continuity of / implies that A q = for all vectors q € Z d . Given an 
irreducible vector q, let bf = ai q for any integer I > 1. Then, the sequence (&/);>i 
is in ^(N) and satisfies 

oo 

vz>i 5>& = o. 

fc=l 

Haar proved that these conditions imply that f>? = for all Z > 1, see [42l Chapter 1, 
no. 129]. This argument holds in all directions q, so the result follows. □ 

We refer to the next section for more general results that explain how to recover 
the coefficients a n from the values A q . 

4. The jump operator 

In order to study the regularity properties of normally convergent Davenport 
series, it is useful to consider the linear operator J which maps the sequence of 
Davenport coefficients to the sequence of jumps, and which is defined by 

V(a„) neZ d € I 1 J ((an)„gz<0 = (A)qez d e (9) 

where the coefficients A q are given by ([7]). The key results concerning this mapping 
follow from those obtained in |35j in the one-dimensional case; this is due to a 
remarkable decomposition that we now present. 

Let I d denote the subset of 7L d formed by the irreducible vectors and let V m de- 
note the vector space of odd sequences (a n ) ng zd that are supported by the multiples 
of such a vector m £ I d , i.e. such that 

supp((a„)„ eZ d) C Zm. 

Any odd sequence indexed by Z d may be decomposed as a sum of sequences b m such 
that b m G V m . However, as the vector subspaces V m and V- m coincide, in order 
to ensure the uniqueness of the decomposition, we shall privilege the irreducible 
vectors whose first nonvanishing coordinate is positive. The set of those vectors is 
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therefore X% = l d n 7L\. As a consequence, we obtain the following unconditional 
Schauder decomposition: 

i x = (K.ru 1 ), (10) 

meaning that any sequence in i 1 may be written in a unique manner as the sum in 
the f 1 sense of an unconditionally summable family indexed by m £ 1+ of sequences 
in V m n t 1 . Moreover, the operator J maps the subspace V m n i 1 to V m n It 
follows that, in order to study J, it suffices to analyze its restriction J m to the 
subspace V m n I 1 , for any fixed vector mElf. 

On top of that, let S m denote the operator of subsampling with step to, which 
is defined by 

S m ({ a n)n& d ) — ( a hn)l>l 

for any odd sequence {a n ) n ^d. As we consider odd sequences only, it is clear that 
the restriction of S m to V m is one-to-one. In fact, we even see that S m maps V m n£ p 
onto £ P {E). 

In the one-dimensional already mentioned above and illustrated by ([3]), 

one assumes that the sequence of Davenport coefficients is supported on N, instead 
of supposing that they form an odd sequence indexed by Z. Thus, the operator 
J has a simpler counterpart which has already been considered in |35j ; this is the 
jump operator J defined by 




J{{b n ) n >i)= \>M) er°(N), (11) 

q>l 

for any sequence (b n ) n >i € ^ 1 (N). The following straightforward lemma shows that 
all the mappings J m essentially reduce to J . 

Lemma 1. Let us consider a vector m G Then, for any sequence a — (a n ) n<£Z d 
in V m l~l i 1 , 

S m (J m (a)) = 2J(S m (a)). 
Therefore, the following diagram is commutative: 

v m n 1 1 v m n £ x 



s„ 



It follows from Lemma [1] that the mapping J can be inverted on each subspace 
V m H £ by means of an inversion formula for the one-dimensional jump operator 
J . This formula has been obtained in 0F| and is recalled in the statement of 
Proposition [5] below. It makes use of the Mobius function [i, which is defined on 
the positive integers by fi(n) — if n is not square-free, and by /x(n) = (— l) fc if 
n is square-free and admits exactly k prime divisors. The inversion formula holds 
on the subspace T(N) of f x (N) that is formed by the sequences (6 n )„>i for which 
the series ~}2, n T(n)\b n \ converges, where r(n) denotes the number of divisors of n; 
this merely means that the restriction of J to that subspace is one-to-one. It is 
well-known that the sequence r(n) grows slower than any positive power of n, in 
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the sense that r(n) = o(n £ ) as n goes to infinity, for all e > 0. This is a plain 
consequence of the fact that 

lim sup l0 f l0gn log rin) = log 2, (12) 

n^oo log n 

see e.g. [3] Theorem 13.12]. This implies in particular that, for any real 7 larger 
than one, T(N) contains the space J" Y (N) of all the sequences b = (6 n )„>i such that 

|&|.F-v(N) = supn 7 |&„| < 00. 

n>l 

Thus, the restriction of J to each ^(N) is one-to-one. The next result even shows 
that J is a bicontinuous automorphism of J 77 (N) . We refer to 35 for its proof. 

Proposition 2. The operator J induces a one-to-one mapping from T(N) into 
^(N). More specifically, for any sequence B = (B q ) q >i in the image set J^(T(N)), 
the equation 

B = J(b) 

admits exactly one solution b — (b n )n>i in T(N), namely, the sequence defined by 

00 

Vn>l b n = J2Kl) B ln- (13) 
1=1 

Moreover, for any real 7 > 1, the operator J induces a bicontinuous automorphism 
of the space J ?7 (N) whose inverse is given by US\) . 

Thanks to Lemma[l] Proposition[2]naturally extends to the multivariate setting. 
In fact, we now establish that the higher-dimensional jump operator J induces a 
bicontinuous automorphism on the space J 77 defined as follows. 

Definition 1. The space J 77 is the vector space composed of the odd sequences 
a = (a„)„ eZ d satisfying 

|a|jr 7 = sup |n| 7 |a„| < 00. 

If 7 is larger than the dimension d of the ambient space, it is clear that J 77 may 
be seen as a vector subspace of I , so that the operator J is well-defined on J 77 . 
In the opposite case, J 77 is not necessarily included in I 1 , However, if 7 > 1, the 
formula ([9]) still has a meaning, because all the series ([7]) converge, which enables 
us to define the operator J on J 77 as well. This may also be seen as a consequence 
of Lemma [TJ along with the fact that 5 m (J 77 ) is contained in ^? 1 (N). 

Proposition 3. For any 7 > 1, the jump operator J is a bicontinuous automor- 
phism of J 77 , and its inverse is given by 



VA = (A q ) qezd G J 77 J- 1 (A) = -V fi(l)A ln . (14) 




Proof. Given 7 > 1, let a = (a n ) nG id be a sequence in J 77 , and let A = {Aq)^^ 
denote its image under J, that is, A = J(a). Then, for each vector q € , 

IA|< 2 | h9 |< 2 |^<^H^, 

where £ is the Rlemann zeta function. Therefore, J is continuous on J 77 with 
operator norm at most 2^(7). 
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In order to study the invertibility of the operator J on J 77 , let us begin by 
observing that, in a way similar to (|10[) . the latter space may be decomposed as 
the direct sum over m £ li of the subspaces V m n J 77 . Moreover, the subsampling 
operator S m induces a one-to-one mapping from V m f~l J 77 , which is included in 
V m n l 1 , onto J-" 7 (N). Then, letting ir m denote the projection onto V m , we deduce 
from Lemma [1] and Proposition [2] that the following diagram is commutative: 



. JT7 



n J 77 ■ 

S m 



2J 



n J 77 
• J" 7 (N) 



This implies the invertibility of the operator J on the space J 77 . In order to obtain 
an explicit formula for its inverse J -1 , let us make use of the diagram to infer that 
the equation A = J(a) implies 

S m (Tt m (a)) = - l 7~ 1 (5 m (7T m (A))) 

for all Tji € where J~ x denotes the inverse of the one-dimensional jump operator 
J on .F 7 (N). By means of (fT3|) . we deduce that 



S m {a) 




»(l)A 



I km 



k>l 



and (fn| follows. Finally, the method that we used above in order to show the 
continuity of the operator J also applies to its inverse because the Mobius function 
is at most one in absolute value; proceeding in this way, we deduce that J -1 is 
continuous with operator norm at most ^("f)/2. □ 

We shall show in Section [6] below that the statement of Proposition [3] may be 
extended to the case where < 7 < 1, up to replacing J 77 by a subspace formed 
by sequences whose support satisfies an additional sparsity assumption. 



5. Pointwise Holder regularity 

The section contains general results on the pointwise regularity of multivariate 
functions with a dense set of discontinuities, a class in which the typical Davenport 
series fall. We begin by recalling the appropriate definitions. 

Definition 2. Let / : M. d — > K be a locally bounded function, xq £ M. d and a > 0. 
The function / belongs to C a (xo) if there exist C > and a polynomial P Xg of 
degree less than a such that for all x in a neighborhood of xo, 

\f(x)-P X0 (x)\<C\x-x \ a . (15) 

The Holder exponent of / at xq is then defined by 



/ l/ (a; ) = sup{a>0|/eC7 Q ( ; r )}. 
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Note that hf takes values in [0, oo]. The following lemma yields an upper bound 
on the pointwise Holder exponent of functions that have a dense set of disconti- 
nuities, and will be applied to Davenport series in the following. It is a direct 
extension to the multivariate setting of Lemma 1 in |33j . 

Lemma 2. Let f : R d — > R be a locally bounded function and let xo £ M. d . Then, 

fa^xo) <liminf 1 1 ° gA/(s) | . (16) 

s-Kc log|s-X | 

Proof. We may obviously assume that hf(xo) is positive, in which case / is con- 
tinuous at xq. Then, given a positive real a less than hf(xo), there exists 6, C > 
and a polynomial P Xo such that (|15[) holds for any x in the open ball B(xo,<5). 

Now, let s be a discontinuity point of /, which thus necessarily differs from xo- 
The magnitude A/(s) of the jump of / at s is positive, as well as e — A/(s)/6. 
Owing to the definition of A/ (s) and the continuity of the polynomial P Xo , there 
exist two points x% and X2 such that 

\f{xi)-f(x 2 )\ > A/(s)-£ and \P Xo (x x ) - P Xo (x 2 )\ < e. 

These two points may be chosen arbitrarily close to s, for instance within range 
\s — xo\/2 from that point. Therefore, at least one of them, denoted by x(s), satisfies 

kW- S |<^ and \f(x(s)) -P X0 ( X ( S ))\>^M. (17) 



Let L denote the right-hand side of (|16|) . which we may obviously assume to 
be finite. Thus, there exists a sequence (s n )n>i of discontinuity points of / which 
realizes the lower limit L; the points s n necessarily differ from Xo but they converge 
to that point. For each integer n > 1, the above procedure yields a point x(s n ) for 
which (fTTl) holds. The resulting sequence (x(s n ))„>i thus satisfies 

,. log|/(^(g»)) -P X0 (x(s n ))\ 

hmsup : t— — r : < L. 

n^oo log \x{S n ) - Xq\ 



In the meantime, (1151) implies that 

liminf Wf(*M) - p *o(*M)\ > a . 

«-+oo log \X{S„) - X \ 

We deduce that a < L, and the result follows from the fact that a may be chosen 
arbitrarily close to the Holder exponent hfixo). □ 

Let us apply Lemma [2] to multivariate Davenport series. For each vector q in 
Z d , let V q denote the set of all integers p £ Z such that gcd(p, q) — 1. We remark 
that a vector q is irreducible, i.e. belongs to X d , if and only if V q contains zero. 
Moreover, the set Hd which indexes the hyperplanes H p q in a unique manner is 
actually formed by the pairs (p,q) such that p £ V q and g e Z|. Now, given a 
point xo G R d , let 5^(xq) denote the distance between xq and the hyperplanes H p ^ q 
with p £ V q , that is, 

s T( x o) = dist [ x , II H p>q I = r^r inf \q-x - p\, (18) 



ist [ x , I) H pq ) = r-r inf 



where |g| is the Euclidean norm of the vector q. It is obvious that 5^ q (xo) coincides 
with <5^(xo), so there is no loss of information in assuming, whenever necessary, that 
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q is in 7A_. Also, it is easy to see that the set V q is invariant under the translations 
of the form p i— > p + k ■ q with k £ Z d , which makes it clear that the function 5^ is 
Z d -periodic. 

The analysis of the discontinuities of Davenport series that we led in Section [3] 
above ensures that every hyperplane H p q indexed by (j>, q) £ %d contains a dense 
set of points at which the Davenport series has a discontinuity of magnitude \A q \, 
with the proviso that the sum A q defined by ([7]) does not vanish. These observations 
then yield the following corollary to Lemma [2J 

Corollary 1. Let f be a Davenport series with a = (a n )„ e zrf £ l x . Then, 



Vzn e M d hf(xn) < liminf 



gGsupp(A) 



log I At I 

log^(^o)' 



where the sequence A = (A q ) qe %d of jump sizes is the image of the sequence a under 
the jump operator J defined by (0). 

In the above statement, we adopt the usual convention according to which the 
lower limit is infinite if the index set supp (A) is finite, in which case the bound is 
trivial. We now illustrate Corollary [1] by pointing out a class of Davenport series 
whose Holder exponent vanishes everywhere, as a direct consequence of the previous 
upper bound. These series are characterized by the fact that the magnitude \A q \ 
of the jumps does not become too small as q goes to infinity along a subsequence 
satisfying particular arithmetical properties. 

In order to specify these properties, let us begin by observing that S^(xq) may 
sometimes be bounded above by l/\q\ infinitely often, up to a logarithmic factor; 
also, in view of the statement of Corollary [TJ we may restrict our attention to the 
vectors q for which A q does not vanish. The situation described above then occurs 
precisely when the support of the sequence A of jump sizes is regular in the sense 
of the next definition, which makes use of the function k defined as follows: For 
every point xq £ R d and every infinite subset Q of Z , 



K 



(xq , Q) — lim sup 



log inf \q-x - p 

logM 



Note that, as a consequence of the periodicity of the function 6^, the function 
k( ■ , Q) is Z d -periodic. 

Definition 3. (1) An infinite subset Q of 7L d is said to be regular if the fol- 
lowing condition holds: 

\/x £R d k{x ,Q)<1. 

(2) A Davenport series with coefficients given by a sequence a £ i 1 is regular 
if supp (J (a)) is a regular subset of Z d . 

It is clear that any infinite subset of a regular set is also regular. Moreover, 
the fact that a set Q is regular roughly means that the sets V q , for q £ Q, do not 
have exceptionally long gaps. In order to elaborate on this remark, let us focus 
on the one-dimensional case and give some heuristic arguments. In that situation, 
the sets V q are gZ-periodic, so that it suffices to analyze their gaps in the interval 
{!,..., q — I}, which clearly amounts to examining the difference between two 
consecutive numbers prime to q. There are <p(q) such numbers, where <f> denotes 
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Euler's totient function. Hence, in the absence of exceptionally long gaps, the 
intervals between two consecutive numbers prime to q would have length of the 
order of q/ <j){q). Thus, the infimum arising in the definition of k(xq, Q) would grow 
at a comparable rate, up to constants, and k{xq, Q) would actually vanish. This is 
due to the fact that q/4>(q) = O(loglogq) as q — > oo; indeed, it is known that 

lim inf < t^Sl. \ g J g q = e _ T j 

9^00 q 

where 7 denotes the Euler-Mascheroni constant, see e.g. [3j Theorem 13.14]. In 
general, though, there may exist exceptional gaps of length much larger or smaller 
than q/4>(q) between the numbers prime to q, so the above arguments are not always 
applicable. The literature seems rather scarce on that difficult topic, apart from a 
series of papers by C. Hooley [251 15CTI [5Tj . 

In some cases, the previous heuristic arguments can be turned into a rigorous 
proof. For instance, let us suppose that Q is the set Z N = {l k , k > 1} of integer 
powers of a given prime number I > 2. For any k > 1, the set V\k is formed by 
the nonmultiples of I, thereby having gaps of length one only. Hence, the infimum 
of \l k xo — p\ over p g is at most two, so that k(xo, l N ) = 0. A one-dimensional 
Davenport series whose coefficients a = (a n ) n >i are supported in Z N will be termed 
as Z-adic. The sequence A = J (a) of jump sizes is then also supported in Z N . As a 
result, any Z-adic Davenport series is regular. 

Note however that in slightly more complicated examples, the above infimum 
may not easily be bounded by a constant, because the sets V q may be chosen in 
such a way that they have longer gaps than above. As an illustration, still in 
dimension one, assume that Q is the set of all primorials of prime numbers, that is, 
the set of all integers qk = Pi ■ ■ ■ Pk for k > 1, where pi is the i-th prime number. 
Then, V qk is the sequence of integers prime that are not a multiple of any of the 
primes pi, . . . ,pk- It follows that V qk has gaps of size at least Pk+i — 1, which tends 
to infinity as k — > 00. 

We now introduce a definition that bears on the asymptotic behavior of the 
sequence indexed by Z d ; we shall apply it in what follows to the sequence formed 
by the jump magnitudes \A q \. 

Definition 4. Let b = {b q ) qe zd be a real- valued sequence and let Q be an infinite 
subset of 7L d . We say that the sequence b has slow decay on Q if 

Uminf - l0g| ^ =Q. 

%q log kl 

Note that this condition is more and more restrictive as the subset Q becomes 
smaller; to be precise, any sequence with slow decay on a given set has slow decay on 
all its supersets. In addition, the above definition is in stark contradiction with that 
of the sets J- 1 that we introduced in Section|3]and on which the jump operator J is 
a bicontinuous automorphism. In fact, it is easy to see that the sequences belonging 
to the sets J 77 , for 7 > 0, do not have slow decay on Z d . 

Combining the previous two definitions, and calling upon Corollary[TJ we readily 
deduce the following result. 

Corollary 2. Let (a n ) n &z d be a sequence in £, and let A — (A q ) qeIi d be its image 
under the jump operator J. Let us assume that A has slow decay on a regular subset 
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ofZ d . Then, 

Vx g R d hf(x ) = 0. 

A typical situation encompassed by the above setting is that of a regular Dav- 
enport series for which the sequence A of jump magnitudes has slow decay. In 
dimension one, this is the case of the Davenport series of the form 

fe=i 

where I is a prime number and a is larger than one. Indeed, being an /-adic 
Davenport series, it is regular. Furthermore, let b = (b n ) n >i denote the sequence 
of its coefficients, namely, b n is equal to k^ 01 if n = l k ° for some integer ko > 1, 
and vanishes otherwise. Then, one easily checks that the sequence of jump sizes 
B = Jib) admits the same expression except that we have to replace k$ a by the 
sum X)fe>fe ^ a m t ne nrs t case. It follows that B has slow decay, and the previous 
result ensures that hf t a (xo) vanishes everywhere. 

6. Sparse Davenport series 

Recall that Corollary [T] above provides an upper bound on the Holder exponent 
of a general Davenport series. In this section, we shall show that, under further 
assumptions, this bound gives the correct value of the Holder exponent, which 
will ultimately enable us to perform the multifractal analysis of the corresponding 
series, see Section [7] Our main assumption implies that there cannot be too many 
nonvanishing terms in the Davenport expansions, and boils down to a sparsity 
condition on the support of the sequence of Davenport coefficients. Note that, in 
one variable, the only case where one can determine the spectrum of singularities 
of Davenport series without additional assumptions on the coefficients is precisely 
the case where one assumes that the frequencies satisfy a lacunarity assumption, 
see [39] . 

6.1. Sparse sets and link with lacunary and Hadamard sequences. For- 
mally, we define the notion of sparse set in the following manner. Recall that 
B(0, R) denotes the open ball of M. d with center zero and radius R. In addition, we 
let # stand for cardinality. 

Definition 5. Let Q be a nonempty subset of E rf . The set Q is said to be sparse if 

Um log#(QnB(0,i?)) =Q 
R->oo log R 

In what follows, we say that an Revalued sequence A = (A n ) n >i is sparse if the 
set of its values, namely, {A„, n > 1} is sparse in the sense of the above definition. 
If there is no redundancy in the sequence, i.e. if it is injective, then the sparsity 
condition suggests that its terms do not accumulate excessively, but rather escape 
to infinity quite fast. Indeed, if A is sparse and injective, it is possible to rearrange 
its terms so as to assume that the sequence (|A„|)„>i is nondecreasing. Then, 
one easily checks that the latter sequence grows faster than any power function at 
infinity, specifically, 

,. log |A„| 

hm — = oo. 

n->oo log n 
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A notable case where the sparsity condition holds is given by the sequences 
(A n ) n >i that are both separated, meaning that 

3C > Vn, m > 1 n^m => |A„ - A m | > C, 

and lacunary, in the sense that 

3C >0 Vn,m>l n^m =► |A„ - A m | > C'(|A„| + |A m |); 

these two notions are standard in the study of nonharmonic Fourier series, see for 
instance [37J El]- Indeed, let us assume that the two above conditions hold, and 
let Afj (A) collect the indices of the terms of the sequence within distance between 
2 J ~ 1 and 2- J from the origin, that is, 

A/}(A) = {n> 1 1 2 j - 1 < |A„| < 2 J '}, 

where j > 1. The lacunarity assumption entails that the open balls B(A„, C"|A n |), 
for n > 1, are disjoint. Therefore, the balls B(A n , C'2^ 1 ) indexed by n € A/"j(A) 
do not intersect either. Meanwhile, all these balls are included in the open ball 
centered at the origin with radius (1 + C j2)2^ . Comparing the volume of these 
balls, we infer that #Aj(A) is bounded above by (1 + 2/C') d . In addition, the set 

M)(A) = {n>l||A n |<l} 

is necessarily hnite; in fact, as a result of the separateness condition, the balls 
B(A„, C/2), for Ao(A), arc disjoint and included in the open ball centered at zero 
with radius 1 + C/2, and the same volume comparison argument implies that 
#7Vo(A) is at most (1 + 2/C) d . As a consequence, the sequence A is sparse. 

The above approach also enables us to write the sparse set {A„, n > 1} as a 
finite union of sets of the form {\^\ n > 1}, where each sequence (X^) n >i is a 
Hadamard sequence, which means that it is separated and satisfies 

\\ {k) I 

3C" > 1 Vn > 1 > C". 

\X ( n k) \ ~ 

Each of these sequences is obtained first by considering the even values of j and 
retaining only one term of the initial sequence A among those indexed by A/} (A), 
and then by handling the odd values of j . Note that a Hadamard sequence is clearly 
lacunary, and therefore sparse, but the converse need not hold. In addition, a finite 
union of Hadamard sequences need not be lacunary. 

6.2. Decay of sequences with sparse support and behavior of the jump 
operator. Recall that the spaces T 1 play an important role in the analysis of the 
jump operator J; in fact, we proved in Section @] that the jump operator J is a 
bicontinuous automorphism of the spaces J 77 , for 7 > 1. We shall now obtain an 
analogous result in the case where < 7 < 1, up to a sparsity assumption. To be 
precise, the set J- 1 will be replaced by the subspace Fg formed by the sequences 
with sparse support, namely, 

where S denotes the vector space of all odd sequences a = (a n ) ne %d for which the 
support supp (a) is a sparse subset of Z d . 

Before studying the behavior of the jump operator on the spaces Fg, let us 
point out that they may be used to characterize the decay of a sequence a € S. 
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Specifically, we measure the rate of decay of such a sequence by considering 

7. = Hminf = *«Kl . (19) 
™^°°, , log n 

TiGsupp(a) O I I 

In particular, a sequence a € S has slow decay on its support in the sense of 
Definition @] if and only if j a vanishes. One then easily checks that 

7 a = sup{7 > | a e J 77 }. 

It will also be useful to remark that, due to the sparsity of the support, "f a may as 
well be seen as a critical exponent for the convergence of a series. Indeed, the fact 
that a is in S implies that 



7a = sup < 7 > 



£ InlKI 1 ^ <oo} =inf<j 7 >0 

nGZ d 



J2 M|0n| 1/7 = 00 



(20) 

As we shall show below, the exponent j a will play a crucial role in the study of 
the multifractal properties of the Davenport series with coefficients given by a, see 
Corollary [3] below as well as the results of Section [7] 

Let us now describe the action of J on the spaces Tg. The next result may be 
seen as a partial extension of Proposition^ to the case where 7 is no more restricted 
to be larger than one. However, it is weaker because it does not discuss invertibility 
properties and the target set of the jump operator J is the intersection 

F 1 -- = p| T 1 ^ (21) 

£>0 

instead of the mere space J 77 . Note that J 77 ' - is endowed with the natural Frechet 
topology inherited from the norms on the spaces J 77-6 . 

Proposition 4. For any 7 > 0, the jump operator J induces a continuous mapping 
in the following way: 

73 — r<- 

Proof. Given 7 > 0, let a — (a„)„ eZ d be a sequence in J 7 ^, and let A — (A q ) q£ %d 
denote its image under J, that is, A ~ J(a). Then, for each vector q G , 

00 00 1 

\AJ<2^\a lq \<2\a\^ {l T^ a)} - 
1=1 1=1 m 

In order to give an upper bound on the last sum, we split the index set into dyadic 
intervals. For each integer j > 0, we have 



2 
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E < ^T#(su P P (a) n B(0, M2^)). 

The support of the sequence a is sparse; thus, for all e > 0, its intersection with 
the open ball centered at the origin with radius |g|2 3+1 has cardinality at most 
C e \q\ e 2 £: ' for some real C e > that depends on neither q nor j. We deduce that 

\A\j*-. = sup \qr e \A g \ < 2|a|^C E E 2(e ~ 7)j = r^M^' 
with the proviso that e < 7. The result follows. □ 
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Note that, even when a sequence a has sparse support, the support of the asso- 
ciated sequence J{a) of jump sizes need not be sparse; this is why the target set in 
the above statement involves .F 7 ~ £ , but not -T 7 ^ -6 • Moreover, the jump operator 
on IFg entails a slight loss in the speed of decay in the sense that the target set is 
not exactly J 77 , as in Proposition [3l but rather J 77 ' - . Still, a simple adaptation of 
the above proof shows that J(a) actually belongs to T 1 when a is a sequence in T 1 
for which supp (a) is composed by the values of a separated and lacunary sequence, 
which is stronger than assuming that a is in S. 

6.3. Pointwise regularity of sparse Davenport series. Now that the notion 
of sequence with sparse support has been defined, we are in position to introduce 
the notion of sparse Davenport series. 

Definition 6. A Davenport series with coefficients given by a sequence a £ I 1 is 
sparse if the support supp (a) is a sparse set, that is, a € S l~l i 1 . 

In order to recover the regularity of the Davenport series / at every point, we 
shall assume, in addition to the sparsity of the support, that there is no cancellation 
in the sums (0 defining the jump operator, in the sense that A q is at least of the 
order of magnitude of its largest term. To be specific, for any sequence a — (a„)„ e ^d 
in i 1 , we may consider the even sequence a = (a q ) qeIl d in £°°(Z d ) given by 

Vq G Z d a q — sup |az g |. 

i>i 

This enables us to define in the following manner a sublinear operator M on I 1 , 
which we refer to as the maximal operator: 

The jump operator and the maximal operator both act on the sequences a = 
(a n ) neZ d in i 1 , and the asymptotic behavior of these actions may be compared 
by means of 

a v log | A, | 

6 a = hmSUp 'TZTTi 

,^oo lOg dg 

where A and a denote the sequences J(a) and M(et), respectively. Incidentally, it 
is useful to remark that 

supp (a) U supp (A) C supp (a) , (22) 

and that supp (A) and supp (a) coincide asymptotically whenever 8 a is finite, in 
the sense that they differ by a finite number of points only. The aforementioned 
assumption may now be expressed by means of the following definition. 

Definition 7. A Davenport series with coefficients given by a sequence a £ i 1 is 
asymptotically jump canceling if 9 a > 1. 

More precisely, assuming that there is no cancellation in the sums defining the 
jumps sizes A q amounts to supposing that 9 a is bounded above by one, i.e. that 
the Davenport series is not jump canceling. The next result shows that, in that 
situation, the upper bound given by Corollary Q] becomes an equality and may 
actually be replaced by an expression that is easier to handle. Specifically, the 
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jump sizes A q arising in the bound may be replaced by the Davenport coefficients 
a n themselves, and the distance 5^(xq) may be replaced by 



\ fcez / 1 1 



S n (x ) = dist ( s , U H k,n ) = o ™f \n ■ Xq — k\, (23) 

which means that we may discard the rather complicated coprimeness condition 
arising in (|Tg|) and discussed in Section [5] 

Recall that the Holder exponent vanishes wherever the Davenport series is not 
continuous, i.e. on the set defined by (jSJ). This set is a union of hyperplanes which 
may be written in the form Dji a ), where for any odd sequence b — (b q ) q £%d, 

D b = (J H Ptq . (24) 

(p,<j)e-Hd 

We may therefore restrict our attention to the points at which the series is contin- 
uous, i.e. outside the set Djr a \. Actually, our approach only enables us to recover 
the Holder exponent of the Davenport series outside the set Dj^ta), which may be 
larger than Djr a \ in view of (|22|) . Yet, this slight restriction will not prevent us 
from performing the multifractal analysis of the Davenport series that are not jump 
canceling, because the two sets then differ by a finite number of hyperplanes only. 
In the next statement, as before, A = (A q ) qeZ d denotes the image of a under the 
jump operator J. 

Theorem 1. Let f be a Davenport series with coefficients given by a sequence a € 
I . Let us assume that the series is sparse and not asymptotically jump canceling, 
that is, 

a G S and 6 a < 1. 
Then, the Holder exponent of f at any fixed point xq € M. d satisfies 

u I \ ^ v ■ r lo SlAI ^ y ■ f lc, gK| 

hf{x ) < limmf H < limmf — — -. (25) 

q eV P 7tA } iog<(x ) n ^;~ M \ogS n (x ) 

Moreover, if x does not belong to D M ( a \, the above quantities coincide and the 
Holder exponent may be computed using either of the following two formulae: 

hf(x ) = limmf H and h f (x ) = lmimf — — -. 

The proof of Theorem Q] is postponed to Section [TU] for the sake of clarity. In 
the course of the proof, we obtain a uniform bound on the Holder exponent, which 
may be seen as a consequence of Theorem [1] and which we state now as a separate 
result for future reference. 

Corollary 3. Let f be a Davenport series with coefficients given by a sequence 
a = { a n) n ei d * n ^ • V the series is sparse and not asymptotically jump canceling, 
then 

Vx eR d h f (x )<j a , 

where "f a is defined by M9\). In particular, the Holder exponent of f vanishes ev- 
erywhere when the sequence a has slow decay. 
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7. Implications for multifractal analysis 

The preceding results will allow us to perform the multifractal analysis of some 
multivariate Davenport series with coefficients a E I 1 . From now on, we assume that 
the series is sparse and not asymptotically jump canceling. We begin by describing 
the size properties of the iso-H61der sets, which are formed of the points where / 
has Holder exponent equal to a given h, specifically, 

E f {h) = {x £ R d | h f {x) = h}, (26) 

for h € [0, oo]. To be precise, we compute the local spectrum of singularities of the 
series /, that is, the mapping 

d f (h, W) = dim H (-B/(/i) n W), (27) 

where W is a nonempty open subset of R d . In the previous formula, dimjj denotes 
Hausdorff dimension, whose definition is recalled in Section 111. 21 The spectrum is 
actually governed by the parameter j a which controls the decay of the sequence a 
and is defined by (|TT)|) . 

In view of Corollary [3l the Davenport series / has Holder exponent at most 
7o everywhere. Thus, all the iso-H61der sets Ef(h), for h > j a , are empty. As 
a consequence, the spectrum of singularities is supported on [0,7 a ], and we may 
restrict our attention to that interval in what follows. In addition, if j a vanishes, 
i.e. when the sequence a has slow decay, then the Holder exponent of / vanishes 
everywhere, so that all the iso- Holder sets are empty, except Ef(0) which is equal 
to the whole space K d . That situation being trivial, we may assume from now on 
that 7 a is positive. 

The analysis below does not cover the case where j a is infinite. Note that this 
case includes that in which the sequence a has finite support. In that situation, the 
Davenport series is a finite sum of piecewise linear functions, thereby being smooth 
except on a locally finite union of hyperplanes where its Holder exponent vanishes. 
If "f a is infinite and the support of a has infinite cardinality, the arguments below 
only imply that the Holder exponent of the Davenport series is infinite Lebesgue- 
almost everywhere in R d , and that the iso-H61der sets associated with finite values 
of the exponent all have Hausdorff dimension at most d — 1. It seems plausible, 
though, that the dimension is exactly d — 1. In what follows, we shall therefore 
assume that j a is finite. 

The local spectrum of singularities of / on the interval [0, j a ] is then given by 
the following statement. 

Theorem 2. Let f be a Davenport series with coefficients given by a sequence a € 
I . Let us assume that the series is sparse and not asymptotically jump canceling, 
and that < 7 a < oo. Then, for any real h G [0,7 a ] and any nonempty open subset 
W ofR d , 

d f (h,W)=d-l + —. 

la 

Note that the spectrum of singularities does not depend on the particular region 
W that is considered, and moreover it is nondegenerate, in the sense that its support 
is not reduced to a single point. Consequently, following the terminology of [36] . 
the Davenport series falls in the category of homogeneous multifractal functions. 

We get comparable results for the singularity sets, which are composed by the 
points where the Davenport series / is continuous and has Holder exponent at most 
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a given h, that is, 

E' f (h) = {x eR d \Dj (a) \h f (x) < h}, 

where Dji a ) is the set of discontinuities of /, given by © and written using the 
notation In addition, we prove that the singularity sets belong to the cate- 

gory of sets with large intersection introduced by K. Falconer [315]. This remarkable 
property essentially asserts that the sets are so omnipresent and large in a measure 
theoretic sense that their size properties are not altered by taking countable inter- 
sections. As a matter of fact, the intersection of countably many sets with large 
intersection with Hausdorff dimension at least a given real s still has dimension 
at least s; this is in stark contradiction with the fact that the codimension of the 
intersection of two subsets is usually expected to be the sum of their codimension, 
as is the case for affine subspaces. Formally, the class of sets with large intersection 
are defined in [26] in the following manner. Recall that a G^-set is a set that may 
be expressed as a countable intersection of open sets. 

Definition 8. For any s £ (0,d\, the class Q s of sets with large intersection with 
dimension at least s is defined as the collection of all G^-subsets E of R d such that 

dim H P| <; n (E) > s 

n>l 

for any sequence (<^i)n>i of similarity transformations of M. d . 

The class Q s is closed under countable intersections and bi-Lipschitz transforma- 
tions, and is the maximal class of G^-sets with Hausdorff dimension at least s that 
satisfies those properties, see [351 Theorem A] for a precise statement. Moreover, 
every set of the class Q s has packing dimension equal to d in every nonempty open 
set, see [26l Theorem D]. In what follows, packing dimension is denoted by dimp; 
we refer for example to |27j for a definition of this notion. 

Restricting to G^-sets will be quite a constraint for us here, so instead of consid- 
ering the classes Q s themselves, we shall work with the extended classes Q defined 
by the following condition: For all E CM. d , 

e e~g s BE' G G s E 1 C E. 

It is clear that the class Q s contains Q s and, in view of [5S1 Theorem C(b)], that 
the G^-sets that belong to Q actually belong to the original class Q s . Moreover, 
the extended class Q naturally inherits from Q s its remarkable properties: Q is 
composed of sets with Hausdorff dimension at least s and packing dimension equal 
to d, and is closed under countable intersections and bi-Lipschitz transformations. 

The next result describes the size and large intersection properties of the singu- 
larity sets of the Davenport series /. 

Theorem 3. Let f be a Davenport series with coefficients given by a sequence 
a = {a n )nei d i n % • Let us assume that the series is sparse and not asymptotically 
jump canceling, and that < "f a < oo. Then, for any real h G (0,7 a ], 

7-1/ / t \ _ -=;d-l+hha 

Ef{h) G y 

and, moreover, for any nonempty open subset W ofM. d , 

dim n (E' f (h)nW) =d-l-\ and dim F (E' f (h) D W) = d. 

la 

We refer to Section [TT] for the proof of the two above theorems. 
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8. Convergence and global regularity of Davenport series 

We will now give a few results concerning the convergence of Davenport series, 
when the sequence of coefficients does not belong to t 1 . In that case, the sum does 
not necessarily belong to L°°, so that Holder pointwise regularity may not be a 
relevant notion anymore. We will mainly consider convergence in Sobolev spaces, 
with both positive and negative indices, which allows us to consider simultaneously 
convergence in spaces of functions or, more generally, distributions. Specific addi- 
tional motivations for this section are supplied in Section [SJ where we show that 
the determination of global Sobolev regularity exponents are a preliminary step 
to either the determination of L q regularity (which is needed for the study of p- 
exponents, see Section l9~4)) or the verification of the multifractal formalism (see the 
beginning of Section [9]). 

8.1. Preliminaries on multivariate arithmetic functions. An arithmetic func- 
tion is traditionally a mapping defined on N and valued in R or sometimes in C. The 
usual multivariate extension deals with functions that are defined on N d , see [T6ll25| . 
In this paper, we consider a slightly different setting, with multivariate arithmetic 
functions defined on Z*. 

A first simple example is supplied by the natural extension to Z d of the divisor 
function already mentioned in Section 2J this extension is still denoted by r for 
simplicity. To be specific, for any to € Z!^, we define t(to) as the number of 
decompositions 

m = In with leN and n G ZJ*. (28) 

It is clear that r(m) coincides with r(gcd(m)), where gcd(m) denotes the greatest 
common divisor of the components of the vector m. With the help of (TT21 . this 
implies that r(m) = o(|to| £ ) as m goes to infinity, for any fixed e > 0. In what 
follows, we shall write indistinctly l\m and n\m when (|28|) holds; with a slight abuse, 
we shall also write / = m/n. 

We will also make use of the extensions to the multivariate setting of other 
arithmetic functions, specifically, the sums of 2-th powers of the divisors. Given 
z e C, recall that the one-dimensional arithmetic function a z is defined by 

a x [m) = ^n 2 , 

n\m 

where the sum bears on the positive divisors of the integer to. In the multivariate 
case, we have to draw a difference between integer and vector divisors. Therefore, 
we define two functions of the vectors to e Z^ by 

a z (m) = \n\ z and a z {m) = 2_,l z - 

It is clear that these two functions coincide on N, and that cto(to) = ao(m) = t(to) 
for all to g Z!^. Moreover, for any zeC, one easily checks that a z (m) = cr z (gcd(m)) 
and a z (m) = \m\ z a^ z (m) = |TO| z er_ 2 (gcd(m)). Given that <J z (l) — l z <7- z (l) for any 
integer I € N, we deduce that 

VmeZf a z (m) = ( -^-] or z (gcd(m)). (29) 
Vgcd(TO)y 
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Finally, recall that the Dirichlet convolution of two arithmetic functions A and 
B defined on N is the arithmetic function A* B given by 

VmeN A*B(m) = £ A(n)B(l). 

(Z,n)£NxN 
ln = m 

Similarly, the convolution of a multivariate arithmetic functions A defined on Z^ and 
a one- dimensional arithmetic function B defined on N is the multivariate arithmetic 
function given by 

VmeZf A*B{m) = £ A{n)B(l). 

ln—rn 

8.2. Davenport expansions vs Fourier expansions. Let us now go back to 
Davenport series. Without any assumption on the sequence a = {a n )n&L d 01 Dav- 
enport coefficients, the right-hand side of ([5]) can be inverted using the multivariate 
Mobius inversion formula, which is an easy extension of the one-dimensional case 
and calls upon the Mobius function \x already used in Section SJ However, we start 
by proving it for the sake of completeness. 

Lemma 3. Let f be a multivariate arithmetic function defined on Zf and let g be 
the multivariate arithmetic function given by 

Vm g Zt g(m) = £ f(n). 

n.\m 

Then, the function f can be recovered from g by f = g * fi, i.e. 

Vm G Zt f(m) = £ g{n)n (^) . 



Proof. For any vector ragZf, we have 

E = O E /(*) = E /(*) E * 

»|m n\m k\n h\m k|n|m 



Let us observe that the integer vectors n£2f satisfying fe|n|m are merely of the 
form n = kl, where I ranges over the divisors of the positive integer m/k. Thus, 
the last sum satisfies 

E^©= E ^(^r) = E A*(o = i{fc=m } - 

"EJj l\(m/k) V / l\(m/k) 

h\n\m 

The last equality follows from the well-known fact that the sum of the Mobius 
function over all positive divisors of a given natural number n vanishes except if 
n = 1, where the sum is equal to one. The result follows. □ 

The next proposition results from applying to ([5]) the above inversion formula, 
and will be useful in the determination of the Sobolev regularity of Davenport series. 
It shows how to recover the Davenport coefficients of a series from the knowledge 
of its Fourier coefficients. 
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Proposition 5. Let f be a Davenport series with coefficients given by a sequence 
a = {a n )nei d * n ^ > an d let { c m) m & d denote the sequence of its Fourier coefficients. 
Then, 

, \ - m ( n \ 

Vn 6 Z + a n = -it > — [i — ) c m . 

Proof. A straightforward consequence of ([5]) is that for all mgZf, 

-Trc m m = ^2 a nn. 

n\rn 

The result now follows from applying Lemma [3] to the arithmetic functions f(n) — 
a n n and g(m) = — ■Kc m m. Note that these functions take values in R d and not 
merely in R. However, Lemma |3] obviously extends to this case. □ 

8.3. Regularity of the sum of a Davenport series. Without any assumption 
on the odd sequence a = (a n ) neZ d of Davenport coefficients, we may define an odd 
sequence (c m ) m &Z d with the help of ([5]). This detour via Fourier series will allow 
us to study the convergence of the Davenport series a n {n ■ x}, even when a is 
no longer assumed to belong to i 1 . Indeed, we shall see that, in many functional 
settings, when the associated Fourier series c m sin(27rm • x) converges, then the 
partial sums of the Davenport series converge to the same limit. 

In order to be more precise, let us begin by recalling that the spaces J- 7 '~ are 
defined in terms of the sequence spaces T 1 by means of (|2Tj) . Moreover, let F 1 
denote the space of distributions whose Fourier coefficients belong to J 77 , and by 
F 1: ~ the space of distributions whose Fourier coefficients belong to T 1 ^ . 

In addition, for any odd sequence a = (a„)„ eZ d, we denote by / the partial 
sums of the corresponding Davenport series, that is, 

f N (x) = ^2 a "{ n ■ x }- 

|n|<N 

The next result discusses the convergence properties of the sequence (/ )n>1 i n the 
spaces F 1 and F 7, ~. A noteworthy consequence lies in the fact that the Davenport 
series a n {n ■ x} converges in the sense of distributions when the coefficients a n 
do not increase faster than any polynomial. 

Proposition 6. Let 7 G R, and let a E J 77 . 

• 1/7 < 0, then the sequence (f N )N>i converges in F 1 '^ to a distribution f 
which belongs to F 1 . 

• ^0 < 7 < 2, then the sequence (f N )N>i is convergent in i^ mln f 1 .7},-_ 

• z/ 7 > 2, then the sequence (f N )N>i converges in F 1- to a distribution f 
which belongs to F 1 . 

Proof. It follows from ([5]) that the Fourier coefficients of the partial sum f N are 
given by 

m 1 \ ^ a, n 



MULTIVARIATE DAVENPORT SERIES 



23 



The condition \n\ < N is necessarily satisfied as soon as N is greater than or equal 
to |m|, so that each sequence (c^)at>i is ultimately constant equal to 

_ 1 \ a n 
° m ~~2^ ^ T' 

(i,n)ez* XZ^ 

ln = m 

Moreover, given that the sequence a belongs to the space J 77 , it is easy to check 
that for all TV > 1 and m € Z», 

which implies in particular that for all m G Z!~, 

|c m | < ^a!_ 7 (m). (30) 
7r|m| 

The following estimates on the one- variable arithmetic functions er z (m) for z S K 
may be found in [45j : 

z<-l (T z (m)= 0(1) 

-1 < z < Ve > cr z (m) = 0(m e ) 

0<z<l Ve > cr z (m) = 0(m z+e ) 

z > 1 c^(m) = 0(m z ). 

Thanks to ([29| . we easily deduce the following estimates on the corresponding 
multivariate functions: 

z<-l cr^m) = 0(1) 

-l<z<0 Ve > cr z (m) = 0(|m| £ ) 

0<z<l Ve > cr z (m) = 0(|m|^ +e ) ( '' 

z > 1 => cF z {m) = 0(\m\ z ). 

It now follows from (|3"TT) that the |c^| satisfy the estimates of Proposition [5] uni- 
formly in N, and their limits \c m \ satisfy the same estimates. 

Convergence in the corresponding function spaces follows immediately by apply- 
ing the same approach to the differences — c m , starting from the observation 
that for any fixed e > 0, 

c m - c m < — g x _ 7+g (m) 
7r|m|iV e 

for all iV > 1 and m G Z^. □ 

Our purpose is now to determine in which Sobolev spaces the Davenport series 
with coefficients in the space J 77 do converge. Let us recall that the Sobolev space 
H s is characterized by the following condition on the Fourier coefficients: A Z d - 
periodic odd distribution / belongs to H s if the sequence (c rra ) meZ d of its Fourier 
coefficients satisfies 

1/1^ = E MVP 5 < 00. 

Note that, if s < 0, this defines a space of distributions. In order to state sharp 
results, we shall also need the following slight modifications of H s . Specifically, let 
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tig be the space of all Z d -periodic odd distributions / whose Fourier coefficients 
satisfy 

\f\h = E i 



< 00, 



and let Hg , and H s 



(1 + log\m\) 5 
be the spaces defined respectively by 



H. 



s,+ 



S+e 



and 



e>0 



e>0 



Before proceeding, let us begin by observing that the Fourier coefficient indexed 

by to = (mi, . . . , T7id) € of the function which maps x = (x\, • • • ,Xd) to {xi} 

is equal to l{ m2= ... =md= o}/ m i- Therefore, this function fails to belong H 1 / 2 but 
l li 

belongs to It follows that, no matter how large 7 is picked, we cannot expect 



1 /2 

substantially better results than convergence in H 1 , . 
less than 1/2, then 



In addition, note that if s is 



|{» 



!}| ff - = ^(2C(2(l-«))) 



1/2 



Thus, for any odd sequence a = (a n ) neZ 
has norm \f\n s bounded above by \ a * 
a consequence, 



Vs < 1/2 



, the Davenport series / defined by @ 
|n| s , up to a multiplicative constant. As 



n < 00 



feH s . 



We will now see how this straightforward result can be improved under the as- 
sumption that a belongs to the space T 1 . We restrict our attention to the situation 
where d > 2, the one-dimensional case being thoroughly studied in 35j. 

Proposition 7. Let us assume that d > 2. Let 7 £ 1, and let a = (a„) neZ ti be a 
sequence in T 1 . Then, the sequence (/ )jy>i converges in the space 



H l,+ 

JJl-d/2- 
H (l+ 7 -d)/2, 



H 



(l+7-d)/2 
1.+ 



ff l/2 



H 



1/2 
1.+ 



if 1 < 

£/ < 7 < 1 

if K 7 < 2 

if 2 < j < d and d > 3 

i/ 7 = d > 3 
ifl>d. 



(32) 



Proof. It follows from Proposition [6] that the sequence (f N )N>i converges to a 
distribution /. Moreover, the Fourier coefficients of / may be bounded with the 
help of pOl). so that 



\f\h < 



IJF2 \ " 
T 2 



|2( S -1) 



^(1+loglTOl)^ 1 -^ 



(33) 
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In order to estimate the above sum, let us split the index set as a union of dyadic 
domains. Specifically, the above sum is equal to the sum over all integers j ' > of 



E 



r( xi- 7 (m) 2 < l — ; ^ m T (y+ 1 )B 7 (y+ 1 ) 



2J <|m|<23 + 1 



(1 + log|m|)' 5 



(34) 



In the previous bound, m 7 and s 7 are defined respectively by 

m 7 (x) = sup cri_ 7 (m) and s y (x) = a\ r (m), 

i-k- | m ,<; 

for any real x > 0. The estimates (|31j) readily imply the following bounds on m 7 : 

'0(x x -T) if 7 <0 
0(x 1 -f +e ) for all e > if < 7 < 1 
0(x £ ) for all e > if 1 < 7 < 2 

0(1) if 7 >2. 
Let us now deal with s-,. For all x > 0, we have 



m 7 (a;) 



*,(*) = £ EH H < E oH M = 

m | <a n »|m |»|<a 

which readily implies that 

{0(^+1-7) if 7 <d 

0(xloga;) if 7 = d 

0(x) if 7 > rf. 

Combining the above bounds, we deduce that the sum in (|34|) is bounded above, 
up to a multiplicative constant, by 

' 2 {A+i{*-i))i / j& if 7 <0 

2 (d+2( s - 7 )+ e )j f or a n £ > if < 7 < 1 

2 (d+2 S -i- 7 + £ )j f or a u e > if 1< 7 < 2 

2 (d+2a-X-7)3 /j<5 
2 {2s-l)j j j5-l 

These estimates are now sufficient to deduce that the limiting distribution / belongs 
to the spaces given by ([32jl. 

On top of that, convergence in the corresponding spaces follows immediately by 
applying the same approach to the differences / — /. As a matter of fact, for any 
fixed 77 > 0, their Fourier coefficients satisfy 

1 n 1 ^ ( 1 + log 1 ?7i 1 y /2 

for all TV > 1 and rnEZf, which implies that (|33|) also holds when replacing |/|fp 
by (l + logATI/^-/! 2 , ■ □ 



if 2 < 7 < d and d > 3 
if 7 = d > 3 
if 7 > d. 
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9. Concluding remarks and open problems 



The study of the local regularity of Davenport series remains a largely open 
field of investigations, with many interesting questions at the crossroad of number 
theory, harmonic analysis and functional analysis; our purpose in this section is to 
list a few of them that we believe of particular interest. A first one consists in the 
verification of the multifractal formalism. We shall not describe this question here, 
because its final and most precise formulation (in terms of wavelet leaders) requires 
the introduction of wavelet methods that go beyond the scope of the present paper; 
we refer to [36 for a mathematical presentation concerning these issues, and to [1] 
for a recent overview on the applications side. The verification of the multifrac- 
tal formalism is a completely open problem for series of compensated pure jump 
functions, whether they be Davenport series (in one or several variables), or Levy 
processes and fields. Let us just mention that Section[5]can be seen as a preliminary 
step in this direction; indeed, a part of this verification involves the determination 
of the Sobolev spaces that contain the function / under consideration. 

9.1. Optimality of Lemma [2j The only cases where we have been able to de- 
termine the exact pointwise Holder regularity of the sum of a Davenport series are 
when the bound given by Lemma [5] is optimal. This is not accidental and actually, 
in all cases of jump functions for which the Holder exponent has been determined, 
it turns out that this bound is optimal: This is the case for Levy processes without 
Brownian component and their extension to the multivariate setting [23l EH 02] , 
for the few cases of Markov processes with nonstationary increments whose multi- 
fractal analysis has been performed [5] , and for the other cases of Davenport series 
which can be worked out [35J [351 I3H] ■ We shall however give below a simple example 
of Davenport series where this is not the case. 

We shall restrict the discussion to the one-dimensional setting, which is easier to 
consider and is sufficient to explain why Lemma [2] is not always sharp, even in the 
setting of Davenport series. Of course, in all generality, (p~6|) is clearly not always 
sharp, as shown by the case in which / is a continuous function, where the bound 
thus obtained is trivial. A natural class of functions for which one might expect 
optimality is supplied by compensated pure jump functions, that is, the functions / 
whose derivative /' in a distributional sense is of the form 



In typical examples, the jump locations x n form a dense subset of the ambient 
space. However, even in the case where / is a compensated pure jump function, 
the bound (|16j) need not be optimal, as shown by the following example of one- 
dimensional Davenport series: 



oo 




(35) 




(36) 



n=l 



where (3 > 3/2. Indeed, the function fp is continuous at zero and jumps at every 
nonvanishing rational p/q written in its irreducible form and the corresponding 
jump has magnitude A/^ (p/q) equal to C(/5) /<Z' 9 ■ Thus, the bound (fl6|) on its Holder 
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exponent at zero is realized by rational numbers of the form 1/q, specifically, 
h h <liminf ° ») , H} = [3. 

We shall prove the following result which shows that this bound is not optimal. 

Proposition 8. Let (3 be a real number larger than 3/2 that is not an integer 
greater than or equal to 3. Then, the value of the Holder exponent of fp at zero is 
given by 

h fp (O)=0-l. 

Proof. Given that the function fp is odd, it is sufficient to study the increment 
fp(x) — fp(0) for positive values of x only. If x £ (0, 1) and n < 1/x, we have 
{nx} — nx — 1/2. Letting |~ • ] denote the ceiling function, we deduce that 

Mz)-Mo) = - m ( x -l)+ 1/ jr 1 ?^< ^ {nx} 



n=l n=[l/x] 



ri/xi-i 



n=l n=\l/x\ 

Let us first assume that j3 < 2. While the first term is merely linear, it is 
easy to see that the second term is equivalent to x l3 ~ 1 /(2 — (3) as x goes to zero. 
Concerning the third term, its absolute value may be bounded by the sum of 1/ nP 
over n > \l/x], which is equivalent to x 13 ^ 1 /(f3 — I). Given that j3 > 3/2, the 
difference \fp(x) — fp{0)\ is thus of the order of and the result follows. 

In the case where /? = 2, the second term in the above decomposition is equivalent 
to xlog(f /x), and the upper bound on the third term is equivalent to x, so the result 
follows as well. 

Finally, let us consider the case in which /3 > 2. The upper bound on the third 
term above is again equivalent to x^~ x / {(i — 1 ). Moreover, the second term may be 
rewritten as 

x ^9=1 = C03 - l)a: - a: J2 ^1' 

n=l n=[l/x] 

and the second term of the latter expression is equivalent to a;' 3 " 1 /(/? — 2). The 
result now follows because j3 is not an integer. □ 

This example opens the possibility of studying the pointwise regularity of Dav- 
enport series for which the bound supplied by Lemma [2] is not optimal. Beyond the 
study of particular functions at particular points, natural general open questions 
are the following. Under simple assumptions on the Davenport coefficients, can one 
show that the bound is optimal at a given point? everywhere? outside a set of 
dimension zero? or almost everywhere? Similar questions can also be raised in the 
more general setting of compensated pure jump functions. 

9.2. Hecke's functions. In dimension one, special attention has been paid to the 
study of very specific Davenport series, namely, Hecke's functions ftp, which depend 
on a parameter j3 € C, and are defined by ([3]). Note that they can actually turn 
out to be distributions when the real part is sufficiently small. 
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These functions have a rich history. They were first considered as functions of 
the complex variable /3, the real number x being merely a parameter. Hecke stud- 
ied their analytic continuation, and his study was later extended by Hardy; these 
results showed that the range of analytic continuation depends on the Diophantine 
approximation properties of the parameter x. As a function of the real variable 
x, the spectrum of singularities was completely determined only in the case where 
3fJ/3 > 2, which leaves open the case where 1 < < 2, see [35] and (|37|) below. 
Note that the counterexample supplied in Section ETTl is closely related with Hecke's 
functions, and we refer to Section 1^731 below for further connections. One could also 
consider multivariate extensions of these functions, specifically, the functions 



where (£ n ) n & d is an °dd sequence taking the values ±1. 

9.3. Spectrum of singularities of compensated pure jumps functions. We 

now go back to the general setting supplied by the compensated pure jump func- 
tions. The examples of such functions whose multifractal properties are known may 
be separated into two large classes. The first class corresponds to the case where 
the jump locations x n appearing in (|35|) are somehow homogeneously distributed; 
this is the case for many examples of Davenport series [351 13H1 13H] or Levy fields 
and processes [33J (Ml IS] ■ The second class is composed of functions for which the 
jump locations x n form a strongly inhomogeneous sequence; such examples have 
been investigated by J. Barral and S. Seuret, and include Levy subordinators in 
multifractal time [9] or heterogeneous sums of Dirac masses [6j [71 [10] . In the het- 
erogeneous case, the obtained spectra strongly differ from those which have been 
exhibited in the present paper. Indeed, they are usually composed of two parts: 
a linear one (for sufficiently small values of the Holder exponent h), followed by 
a strictly concave one. Note however that spectra of a different kind have been 
obtained in [5j; for some Markov processes which differ from Levy processes, one 
meets spectra that are a superposition of linear functions with different slopes. 

In the homogeneous case, the spectra that have been met up to now are lin- 
ear. However, this is not a general rule, even in the particular case of one-variable 
Davenport series. Let us consider for instance the function fp defined by ([3rJ| . 
and suppose that f3 is a noninteger real number larger than two. The local spec- 
trum of singularities of the corresponding Hecke function f)p defined by ([3]) is then 
supported by the interval [0,/3/2] and satisfies 



this follows readily from the approach employed in [35] in order to compute the 
global spectrum of Sjp, which corresponds to the case where W is equal to the 
whole real line. As shown by Proposition^ subtracting the term (([3){x} to Hecke's 
function ij p (x) shifts the value of the Holder exponent at the integers from zero to 
13 — 1. As a consequence, the local spectrum of singularities of the resulting function 
fp is now supported in the set [0,/3/2] U {/3 — 1}. Moreover, the function fp still 
satisfies ([37)1 but, rather than being empty, its iso- Holder set Ef (0 — 1) is equal to 
Z. Therefore, for any open subset W of M that contains an integer, 




V/ie[0,/3/2] WW t^0 open 



d^(h,W) = — ; 



(37) 



d h (0-l,W) = Q. 
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In particular, the local spectrum of fp depends on the particular region that is 
considered. Thus, unlike the corresponding Hecke function S)p, the function fp 
is not a homogeneous multifractal function. Furthermore, unlike that of ftp, the 
global spectrum of singularities of fp is not a linear function. More specifically, 
the graph of this spectrum is the union of a segment and an isolated point. The 
multifractal properties of fp may be put in comparison with those of Riemann's 
function ^ n sin(7rn 2 a;)/n 2 whose global spectrum of singularities has exactly the 
same shape and which is a homogeneous multifractal function, see |32j . 

This example raises several questions. Is there a simple condition on the coeffi- 
cients of a Davenport series which ensures that its spectrum is linear? What is the 
general form of the spectrum of a Davenport series? Similar questions can also be 
asked in the more general setting of compensated pure jumps functions. Note that 
some results on these problems have been obtained by J. Barral and S. Seuret in 
the slightly different setting supplied by the large deviation spectrum, see [5]. 

9.4. p-exponent. Holder pointwise regularity is defined only for locally bounded 
functions, which explains why we always made the assumption that the sequence 
of Davenport coefficients belongs to £ when studying Holder regularity. However, 
Proposition [7] shows that, even when the sequence of Davenport coefficients does 
not belong to l l , and therefore convergence in L°° is no more guaranteed, one can 
obtain convergence in L 2 (this corresponds to the Sobolev space H s considered in 
Section [5] in the case where s is zero), and also for larger values of p; indeed when 
s is positive, the Sobolev embeddings imply that, if / belongs to H s , then it also 
belongs to LP for all p smaller than the critical value po defined by the condition 

1 _ 1 s 
Po 2 d' 

Note also that the specific case of L 2 convergence of one-variable Davenport series 
has already been considered, see jT5] [35] ■ 

In such situations, one can still perform a pointwise analysis of regularity, by 
using a definition of pointwise smoothness which is weaker than Holder regularity 
and is compatible with functions that are not locally bounded: It is the notion 
of TP(xo) regularity, which was introduced by Calderon and Zygmund in 1961, 
see [18] . The next definition is an adaptation of Definition [2] to that setting. 

Definition 9. Let / be a tempered distribution on R d , let p € [1, oo), let a > —d/p 
and let xq G R d . The distribution / belongs to T p (xq) if it coincides with an LP 
function in the open ball T5(xq,R) for some real R > 0, and if there exist a real 
C > and a polynomial P X() of degree less than a such that for all r <E (0, R], 

i\l \f(x)-P X0 (x)\Pdx) <Cr a . 

The p-exponent of / at xq is then defined as 

h*(x ) = sup{a > -dip | / e TP(x )}. 

Note that the Holder exponent corresponds to the case where p = oo, and the 
condition on the degree of P Xo implies its uniqueness. This definition is a natural 
substitute for pointwise Holder regularity when functions in Lf oc are considered. In 
particular, the p-exponent can take values down to —d/p, thereby allowing to take 
into account behaviors which are locally of the form l/\x — xo\ a for a < d/p. 
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Furthermore, similarly to (|26|) and (|27|). we may define the analogs of the iso- 
Holder sets and the local spectrum of singularities by 

E p f (h) = {x e K d | h p f (x) = h} and d p f (h, W) = dim n (E p f (h) n W), 

the latter quantity being referred to as the local p-spectrum of the distribution /. 

Both in the univariate and the multivariate case, the subject of determining 
the p-exponents and the p-spectrum of a Davenport series with coefficients not 
belonging to I 1 is completely open. 

9.5. Directional regularity. The notion of Holder pointwise regularity given in 
Definition [5] does not take into account directional regularity but yields the worst 
possible regularity in all directions. Therefore, all the results obtained in the present 
paper do not take into account possible directional irregularity phenomena. 

We now briefly discuss the notion of directional regularity. Let / be a locally 
bounded function defined on M. d . In order to take into account directional behaviors, 
it is natural to define the Holder regularity at xq in a direction u € R d \ {0} as 
the Holder regularity at zero of the univariate function t i-> f(xo + tu). This 
definition has several drawbacks which stem from the fact that the latter function 
is defined as the trace of / on a line, which is a set of measure zero, see [37] for a 
detailed discussion. Let us now give the definition of anisotropic smoothness which 
is currently used, see e.g. [TT1I57] . 

Definition 10. Let / be a real-valued function defined on R d and bounded in a 
neighborhood of a point xq S R d . Let e = (e%, . . . , ed) be an orthonormal basis of 
M. d and let a — (ai, . . . , ay) be a <i-tuple of nonnegative real numbers such that 
ot\ > ... > ad- The function / belongs to C a (xo, e) if there exist a real C > and 
a polynomial P XQ such that for all x in a neighborhood of xq, 

d 

\f(x)-P Xo (x)\<Cj2\(x-xo)-ei\ a *. 

i=l 

Since, by construction, Davenport series display jumps along hyperplanes, thereby 
being extremely anisotropic by nature, a natural question is to determine their 
pointwise anisotropic regularity; the same remark is also relevant for Levy fields, 
which present the same type of anisotropy, see |24j . These two examples would 
certainly be natural candidates to test possible definitions of anisotropic spectra of 
singularities. 

Finally, note that an extension of pointwise smoothness combining anisotropy 
and the T^(xo) condition is proposed in |37j . This notion could be relevant in 
order to perform the study of the anisotropy of multivariate Davenport series with 
coefficients not belonging to t l . 

10. Proof of Theorem Q] 

Let us begin by comparing the two lower limits appearing in the statement of 
the theorem. First, given that 9 a is bounded above by one, we have 



liminf -= , ; < liminf 



, eBUP p(A) log^(xo) qe lZ7^ log^(zo)' 

where A = J(a) and a = M(a). In addition, replacing 8^ (xq) by S q (xo) in the 
right-hand side does not change the value of the lower limit. Indeed, for any point 
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q in the support of the sequence a, the distance S g (xo) is reached on a hyperplane 
of the form H pq with p £ Z, so that 

x , s |g • ^0 ~ Pi W ■ X Q-P'\ ^ c7V s 

kl kl 

Here, p' — p/r and = q/r, where r denotes the greatest common divisor of the 
integer p and the components of the vector q, so that p' G V q > . Since the multiples 
of q are also multiples of q' , the supremum over all integers I > 1 of \ai q \ is at most 
that of \ai q i\. As a consequence, for all q G supp (a) large enough, there exists a 
point </|g such that 

logy < loga 9 < logQg 
log<5£(x ) ~ log(5 (a;o) ~ logJ^xo)' 

so the lower limit featuring S^xo) in the denominator coincides with that featuring 
S q (xo). Furthermore, \a q \ is obviously bounded above by a q , so that 

log a? < log kg I 
\ogS q (x ) ~ logd q (x )' 

The above discussion, combined with Corollary[TJ finally leads to (f^5|) . In particular, 
since 5 n (xo) is bounded above by l/|n| regardless of the value of Xo, we deduce the 
next uniform bound on the Holder exponent: 

hf(x ) < 7a, 

where j a is defined by (fTi?)) . The remainder of the theorem then follows when the 
sequence of Davenport coefficients has slow decay, i.e. when j a vanishes. 

In order to finish the proof of the theorem, we thus may assume that "f a is 
positive. It remains us to establish that 

h f (x ) > hminf Z 08 '"" 1 . (38) 

for any point xq that does not belong to the set Z?M(a) ■ Lrt us consider an integer 
jo > and a point x eR d such that 2 _ ( j0+1 ) < \x — x \ < 2~ Ja . Since the sequence 
a is in I , we may define 

^x a ,x{Z) = ^ «n {{n ■ X} - {n ■ Xq}) 

nez 

for any subset Z of Z d . Since the Davenport series converges normally, it is clear 
that its increment between xq and x may be written in the form 

fix) - /(x ) = E X0 , x (Z d ). 

Therefore, it suffices to handle the series E xo iX (Z) for Z ranging over a collection 
of sets that form a partition of 7L d . 

To be specific, let us begin by giving an upper bound on |£a; 0j:E (Cj (1 Z d )|, where 
Cj is the domain defined by 

Cj = {xe R d | 2 J < \x\ < 2 1+1 } , 

for any integer j > 0. For every fixed 7 € (0,7 a ), there exists a constant C 7 > 
such that |a„| < C 7 |n|" 7 for all n £ Jr. For the sake of simplicity, we merely write 
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| On | "C |n| 7 in such a situation, thereby making use of the Vinogradov symbol. 
Accordingly, 

{Cj n z d )\ < E Kl < {nfEsupp(a)} ■ 

Furthermore, for every fixed e > 0, we deduce from the sparsity assumption bearing 
on the sequence a that #(supp (a) n Cj) <C 2 £J . As a consequence, 

E l^r 71 {n G su P p(a)} < 2-^'#(supp(a) nCj) « 2(- 7 +^. 

n6CjnZ d 

The union over all integers j > jo is equal to the complement in R d of the open 
ball centered at the origin with radius V a . Thus, summing over all these values of 
j, we obtain 

|£ Xo , x (Z d \B(0,2*>))| «2(- 7 +^° « \x-xor E . 

Now, let J\f xo ,x denote the set of all points n E Z d for which there exists an 
integer k G Z satisfying either n-xo<k<n-xoTn-x<k<n-xo, meaning 
that the hyperplane Hk t n separates the points xq and x, and let Af£ x denote its 
complement in Z d . For each n 6 7VJ„ xl it is clear that n ■ xo and n ■ x have the 
same integer part, so that 

S v (^ i nB(0,2'»)) = ( I -io)' £ a n n. (39) 

n|<230 

If 7 a is smaller than or equal to one, the modulus of the sum in (|39p may be bounded 
above by the sum over j € {0, . . . , jo — 1} of 

£ KH < E M 1-71L {«6 S upp(a)} 

nGCjnZ d neCjDZ d (40) 

« 2< 1 -^''#(supp (a) nCj) « 2< 1 - 7+£ W, 

which entails that 

|£* o ,*(A/^ x nB(0,2->°))| « Ix-xol^ 1 -^^ < \x-xoV- £ . 

In the opposite case where 7 a is larger than one, we assume that 1 < 7 < 7 a and 
rewrite (|39p as a difference of two terms in the following form: 

s *o,*(^ ,x nB ( > 2J ' )) = (x-xa) ■ E a„n-(x-x )- E a n n. 

|n|>2J0 

The first series is normally convergent because, in view of (|40j) . we have 

00 00 

E Kn|<E £ \a n n\^Y, 2(1 ^ +E)3 <°°- 

neAA; 0iX i=0n£C 3 nZ d J=0 

In order to handle the second term, we make use of (|40|) again; this implies that 

00 00 

E < E E i a « n i « E ^ i+e)] « ii - x r i+7 - £ . (41) 

|7l|>2« 



MULTIVARIATE DAVENPORT SERIES 



33 



We finally deduce that 



<€. \x — x 



7— e 



It remains us to consider the behavior of ^x ,x on the set M Xo ,x H B(0, 2- 70 ). Let 
a denote the lower limit in the right-hand side of ()38p . which we may assume to 
be positive, and let a' € (0, a). By definition of a, we have 8 n (xo) > leinl 1 /" for 
n G Z d sufficiently far from the origin. In addition, a n necessarily vanishes when 
the distance 6 n {xo) is zero. As a matter of fact, in that situation, xq belongs to 
a hyperplane -Hfc jTt with k G Z. This hyperplane is represented by a unique pair 
(p, (?) G Hrf, and then p and g divide k and n, respectively. However, xq does 
not belong to I?M(a)j so the index g cannot belong to the support of the sequence 
M(a) = a, which entails that \a n \ < a q — 0. The upshot is that there exists 
a real C Q < > such that 5 n (x Q ) > Ca^a^ 1 / for all n G Z d , which we write 
I On | "C (5 n (a;o) a still using the Vinogradov symbol. As a consequence, 

ISxo^CA/'xoz nB(0, 2 3o ))\ < ^2 Kl < {nEsupp(a) } • 



<2J0 



|r.|<2* 



Moreover, when n belongs to M XOtX , we have \n ■ Xo — k\ < \n ■ (x — xq)\ for some 
integer k G Z, so that 5 n (a;o) is bounded above by \x — xq\. Hence, 

Sn(xo) a 't {nesuMa}} < |x-x o r'#(supp(a)nB(0,2 J «)) « \x - x \ a '2 ejo , 



where the last bound follows from the sparsity assumption bearing on the sequence 
a. We deduce that 

|E*o,»(JV*o,*nB(0,2*'))| « \x-x \ a '~ £ . 
The above approach also enables us to write that 

— , , — , , jo ~ 1 

\a n n\ < \x-x \ a> Y ^^{nesuppCo)} < \x~x \ a ' ^ 2:, #(supp(a) nC 3 ), 

3=0 



|n|<2J0 



n|<2J0 



where the last sum is bounded by 2^ 1+E ^° up to a constant, in view of the sparsity 
of the sequence a. In addition, when 1 < 7 < j a , the bound given by pij) still 
holds when A/£ x is replaced by Af Xo , x - It follows that 



which readily implies that 



-1+7-e 



+ \X - X \ 



-l+a'-s 



Sx 0l x(^o,xnB(Q > 2»'°))-(a;-xo)- ^ 



<€. \x - x 



7— s 



+ F - xq| 



Combining all the previously obtained bounds, we finally get 
|/0r) - f(x )\ « |z - x | 7 ~ £ + \x- xo\ a '~ e 
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when 7 a is smaller than or equal to one, and 

f(x) - f(xo) - (x - xq) ■ 2J a n n < I a: - x \ J ~ E + \x- x a \ a ~ e 

when 7 a is larger than one. In both cases, it appears that the Holder exponent at 
Xq of the Davenport series / is at least the minimum between 7 — e and a' — e. The 
bound ([38| finally follows from letting e go to zero, 7 to 7 a and a' to a. 

11. Proof of Theorems E] and [3] 

Throughout the section, / denotes a Davenport series with coefficients given 
by a sequence a = (a n ) n£ %d in (r. We assume that the series is sparse and not 
asymptotically jump canceling, and that j a is both positive and finite. 

11.1. Locations of the singularities. The first step to the proof of Theorems [2] 
and [3] consists in observing that the iso-H61der sets and the singularity sets of the 
Davenport series / may be expressed in terms of the sets L a (a) of all points that 
are at a distance less than la,^ 1 /" from a hyperplane Hk, n defined as in ((6J for 
infinitely many points n in the support of the sequence a. Put another way, a point 
x belongs to L a (a) if and only if the distance 5 n (x) defined by (|2"3")l is less than 
Icinl 1 /" infinitely often. To be more specific, for any real a > 0, the set L a (a) is 
defined by 

L a (a) = {xeR d \\n-x-k\<\n\ \a n \ 1/a for i.m. (k,n) eZx Z d }, (42) 

where i.m. stands for "infinitely many" . It is easy and useful to remark that the 
mapping a > L a (a) is nondecreasing. 

The connexion between the iso-H61der and singularity sets, and the sets L a (a) 
is now given by the next result. It is a direct consequence of Theorem [TJ along 
with the discussion made in Section [3] above according to which the Davenport 
series / is discontinuous on Djt a \, thus having Holder exponent zero thereon. In 
its statement, M(a) denotes the image of the sequence a under the action of the 
maximal operator. 

Lemma 4. Let h 6 [0,7 a ]. Then, 

E f (h) \ D J(a) C E' f (h) C (D M{a) \ D J(a) ) U p| L a (a). (43) 

a>h 

Moreover, Ef(0) D -Dj(a) and for the positive values of h, 

E' f (h) D L a (h) \ A/(a) and E f (h) D E' f (h) \ (J L a {a). (44) 

a<h 

Lemma[l]suggests that the proof of Theorems [5] and [3] will follow from a detailed 
understanding of the size and large intersection properties of the sets L a (a). This 
is the purpose of the next subsection, but let us just point out here that 

P L a (a) =R d . (45) 

Ct>7a 

Indeed, it is plain that S n (x) < l/\n\ for every n € Z d and every point x € M d , 
and that |a„| 1/,a |n| > 1 infinitely often, when a > j a . We may thus restrict our 
attention to the case where a < "f a in the study of the size and large intersection 
properties of L a (a). 
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11.2. Size and large intersection properties of the sets L a (a), connection 
with the Duffin-Schaeffer and Catlin conjectures. We shall investigate the 
size properties of the sets L a (a) by estimating their Hausdorff measures for specific 
gauge functions. We call a gauge function any continuous nondecreasing function 
g which is defined on [0, e] for some e > and vanishes at zero. The Hausdorff 
measure associated with such a gauge function is then defined by 

oo 

W{E) = hm t H$(E) with U% (25) = ^nf^ j>(M), 

for any subset E of R d . Here, the infimum is taken over all sequences (£/j)i>i of 
subsets of R d satisfying E C [Zj and |£/j| < 5 for all i, were | • | denotes diameter. 
It is well-known that H 9 is a Borel measure on R d , see e.g. [44]. Moreover, the 
Hausdorff measure associated with the gauge function r 4 r s is called the s- 
dimensional Hausdorff measure and is denoted by recall that such measures 
enable one to define the Hausdorff dimension of a nonempty set E C M. d by 

dim H E = sup{s e (0, d) \ H S (E) = oo} = inf{s e (0, d) | H S (E) = 0}, 

see Falconer's book [27] for instance. 

Theorems [2] and [3] state that the iso- Holder and the singularity sets of the Dav- 
enport series / all have Hausdorff dimension between d — 1 and d. Therefore, on 
our way to the proof of these results, we may restrict our attention to the gauge 
functions of the form r i->- r d-i+s^ w ^h < s < 1, as well as slight corrections 
thereof. These corrections are obtained by replacing r s in the previous expression 
by more general functions, specifically, the continuous nondecreasing functions ip 
defined on [0,e] for some e > which vanish at zero, for which r i— > ip(r)/r is 
nonincreasing and positive, and for which the limit 

log^(r) 

s v = hm — 

r->o logr 

exists (this limit is then between zero and one). The collection of all such functions 
is denoted by and clearly contains the functions ri->r s , for < s < 1. For any 
ip £ <fr, it is now plain that the function r M> r d ^ 1 (p(r) is a gauge; the corresponding 
Hausdorff measure is denoted by 'H d ~ 1 ' ip , and the value that it assigns to the set 
L a (a) is discussed in the next statement. 

Lemma 5. For any real a > and any function ip G $, 

J2 M^(KI 1/Q ) <oo n d -^(L a (a)) = o. 

Proof. Let p n = lan] 1 ^ 01 for any n € Z d , and let us consider two real numbers 
A > 1 and 8 G (0, 1]. Let us assume that the series appearing in the statement 
of the lemma converges. So, there necessarily exists an integer T)o > 1 such that 
4p„ < 5 for all n € Z d with \n\ > 770. Then, for any 771 > 770, 

L a (a) n B(0,A- 1) C |J |J {xeB(0,A) | distfo < p„} . 

nesupp(a) fc£Z 
|tl|>f/l |fc|<A|7i 
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Moreover, each set in the union above may be covered by (212 AVd/ p n \) d 1 open 
balls with radius 2p n . Therefore, 

C"" (r) (M«)nB(0,i-l)) < J2 2A\n\(^) (4p„) d -V(4p„) 

-E-nppW V PU ) 

M>*!1 

< SAilGAVdf- 1 WM\an\ 1/a ). 

M>*ji 

Letting 771 00 and 5 -» 0, we deduce that H d ~ 1 ' ,p (L a (a) H B(0, A- 1)) = 0. This 
holds for all integers A > 1, so the result follows. □ 

In particular, letting tp be the identity function in the statement of Lemma [5] 
and letting C d be the Lebesgue measure in R d , we deduce that 

MK| 1/Q <oo £ d (I n (a)) = 0. (46) 

Owing to the alternate expression (|20|) of 7 Q , it is easy to see that the above series 
converges once a is less than j a . Owing to Lemma [4] and the fact that -Dm(o) is a 
countable union of hyperplanes, we deduce that the iso- Holder sets Ef(h) and the 
singularity sets E'^(h) have Lebesgue measure zero when h < j a . Together with 
Corollary [31 this implies that the sets Ef{^ a ) and E'^{^ a ) both have full Lebesgue 
measure. Therefore, hf(xo) = 7 a for £ d -almost every Xo € K d . 

To proceed with the proof of Theorems [5] and El we shall need a kind of converse 
to Lemma[5l which ensures that L a (a) has a positive Hausdorff measure for specific 
gauge functions. The study of the size properties of various classical sets arising 
in the metric theory of Diophantine approximation suggests that such a converse 
should look like the following: 

\n\tp{\a n \ 1 / a )=oo VVKopen V.^ 1 ^ (L a (a) n W) — 'H d ^ 1 ' v (W), (47) 

see for instance [THIH] and references therein. Given that L a (a) is of the form 

K d (ip) = {x G R d I I ra • x - k\ < ip(n) for i.m. (k,n) eZxZ d }, 

where ip : Z d — > [0, 00) is a multivariate approximating function, and thanks to the 
mass transference principle of [12] and the slicing technique of [13] , this would follow 
from the next general statement: The set Kd(ip) has full Lebesgue measure in K d if 
the series ip(n) diverges. However, such a statement is known to be false and we 
refer to Section 5 in [14] for a counterexample. The expected result is actually given 
by a generalization of the Catlin conjecture to the case of dual approximation, which 
has been formulated by V. Beresnevich, V. Bernik, M. Dodson and S. Velani [14], 
and consists in replacing the previous series by 



where | • |oo denotes the supremum norm and <j>d(o) is the number of positive integers 
less than or equal to |g|oo which are coprime with the components of q. We refer 
to [14] for a motivation of this conjecture, and for its relationship with the dual 



MULTIVARIATE DAVENPORT SERIES 



37 



form of the famous Duffin-Schaeffer conjecture. The upshot is that it seems rather 
difficult to provide a converse to Lemma [5] in the form (|47l) . 

As shown by the statement of Theorems and [3J we ultimately describe the 
size of the iso-H61der and singularity sets in terms of Hausdorff dimension, rather 
than using general Hausdorff measures. Thus, we do not need to call upon such 
precise results as those mentioned just above. In fact, we only need to prove that 
appropriate modifications of the set 

L a ( 7a ) = [ x £ R d | 1 7i • x - k\ < \n\ \a n \ lha for i.m. (k,n) eZx Z d }, 

obtained by letting a = -f a in PS)) , have full Lebesgue measure in M. d . Note that 
the set L a (a) has full Lebesgue measure in M. d when a > j a due to (|45)) . and 
Lebesgue measure zero when a < "f a by virtue of Lemma [5] The study of its 
Lebesgue measure for the critical value a — -f a is more delicate. Indeed, if true, (|T7|) 
would imply that L a {^ a ) has full Lebesgue measure when the series ^ n lanl 1 / 7 " 
diverges. However, the coefficients of the Davenport / series may be chosen in 
such a way that the series converges, e.g. when the nonvanishing coefficients are 
given by \a\ m \ — (m 2 |A m |)~ 7 for some positive real 7 and some sparse injective 
sequence (A m ) m >i, in which case L a (j a ) has Lebesgue measure zero, by Lemma[5] 
This means that we shall have to slightly reshape this set in order to ensure that 
we work with a set with full Lebesgue measure. Actually, as shown by the next 
lemma, a slight modification of -E/ a (7a) enables one to recover the whole space. This 
modification is written in the form 

4 P '%a) = {xeM. d \\n-x-k\< H^KI 1 /^) for i.m. (k,n) E Zx A/i}, 

where ip and i are appropriately chosen in <E> and {1, . . . , d} respectively. Here, Aft 
denotes the set of all points n — [n\, . . . , nj) in ^ d such that |n|oo = |7ij|; note that 
the sets A/i obviously form a covering of Z d . Whereas the function tp is crucial in 
order to enlarge the set L a (^ a ) and then to recover the whole space, the index i, 
which is used to retain only some specific frequencies among the support of a, is 
introduced merely for technical reasons appearing in the proof of Lemma [7] below. 
In the following, denotes the collection of all functions <p £ $ with = 1 for 
which at least one of the sets La?' 1 (ja), ■ ■ ■ , £0 (7a) has full Lebesgue measure 
in K d . The next lemma shows that is nonempty. 

Lemma 6. There exist an index i+ and a function (p* with s Vt — 1 such that 

\n\ Lp+(\a n \ 1 / la ) > 1 for i.m. 71 SA/"^. 

In particular, the set La V *'**\"/ a ) is equal to the whole space , and tp+ is in 

Proof. Let (A m ) m >i denote an enumeration of the support of the sequence a. Given 
that the Davenport series is sparse, the sequence (A m ) TO >i is sparse and injective 
and, up to rearranging its terms, we may assume that the sequence (|A m |) m >i 
is nondecreasing. Then, let p m — \a\ m \ 1 ^ la and u m = l/|A m | for any integer 
m > 1. The case in which the sequence (u m / p m )m>i does not diverge to infinity 
is elementary. Indeed, in that situation, there exists a constant C > such that 
u m < Cp m for all m belonging to some infinite subset M. of N. Then, the function 
defined by (p*(r) = Cr clearly satisfies the required properties, and it suffices to 
choose i+ in such a way that A/i t contains infinitely many points A m with mcM. 

We may therefore suppose from now on that (u m / p m ) m >i diverges to infinity. 
The definition (fT9j) of j a , combined with the observation that Z d is covered by the 
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sets Ni, ensures the existence of an index such that 

logu m 
lim sup = 1 . 

CZZ. logPm 

In addition, both u m and p m tend to zero as m — > oo. Thus, we may find a sequence 
of indices {m,k)k>i in N along which all the following properties hold: The sequence 
(u mk / Pm k )k>\ diverges to infinity monotonically; log u mk / log p mk tends to one as 
k goes to infinity; for all fc > 1, 

Pm k + 1 < (Pm k ) 
u rn k + l 5: i u m k ) ■ 

It is now straightforward to check that any logarithmic interpolation of the points 
(p mk , u mk ) yields a suitable function ip*. To be specific, any function ip* defined on 
[0, oo) for which 

log^(r) = logu mk - — — — (logw mfc - logu mfc+1 ), 

iOg Pm k iOg Pm k+1 

for all r £ (p mk+1 , Pm k ] and k > 1, clearly belongs to the set $ and meets all our 
requirements. □ 

For any a £ (0, j a ], we now define a mapping T a on the set by letting 

T Q v? : r h> ^(r a/7 ») 

for any function tp £ All the functions r a (^ belong to $ and satisfy st qV = 
ct/"f a , so they roughly behave like r a l la near the origin. Furthermore, T 7a is 
the identity mapping. Now, recall that Lemma |6] yields a function 93* for which 
\n\ (^(letnl 1 / 70 ) > 1 infinitely often; the series ~J2 n \n\T a (p ± {\a n \ i ^ c ') thus diverges. 
We are in a situation where the assumption of Lemma [5] fails and, in fact, the set 
L a (a) does not necessarily have a vanishing T-L d ~ 1 ' Tc ' v *-ma,ss. On the contrary, it 
should be regarded as large and omnipresent in K d in terms of (d — 1 + a/j a )- 
dimensional Hausdorff measure, in the sense that it belongs to Falconer's class 
gd-i+a/f a Q f se ^ g w j^h large intersection defined above. 

This follows from our next result, namely, Lemma [71 which describes the large 
intersection properties of the sets L a (a). The properties are expressed by means of 
the classes G 9 (W) that were introduced in [21] in order to extend Falconer's classes 
to general gauge functions g and open sets W C R d , with a view to establishing a 
suitable framework to describe precisely the large intersection properties of various 
sets arising in the metric theory of Diophantine approximation. In what follows, we 
shall restrict our attention to the classes G d_1 ' ¥3 (R <i ) of sets with large intersection 
in the whole space M. d with respect to gauge functions of the form r i— > r d ~ 1 </?(r) 
with We refer to [21] for a precise definition of those classes and a description 

of their main properties, and we content ourselves here with recalling that the class 
G d_1, ' /> (]R d ) is closed under countable intersections and bi-Lipschitz mappings, and 
is formed of G^-subsets E of W l satisfying 

VW open H d ~ 1 ^(E n W) = 00 (48) 

for any function tp £ $ growing faster than tp at zero, in the sense that ip/<p 
tends to infinity monotonically, in which case we write tp -< (p. A straightforward 
consequence of these properties is the fact that, when d — 1 + s v is positive, the 
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class G d ~ 1,v (M. d ) is included in Falconer's class g d - 1 + s <c . Let us now describe the 
large intersection properties of the sets L a (a); in the next statement, $ a denotes 
the collection of functions satisfying tp -< T a ip± for some € <!>*. 

Lemma 7. for cmy rea/ a € (0, 7J, 

L a (a) £ p| G d_1 ' ¥ '(M' i ) C ^-i+«/7a_ 

Proof. Let us consider a function ^ G $ Q , and a function <y9* G <£>* for which 
</? -< T a ip+. We shall make use of a ubiquity result, which enables one to deduce the 
large intersection properties of the set L a (a) from the sole fact that a corresponding 
enlarged set, namely, one of the sets L a v (7 ), • ■ ■ , La (j a ) has full Lebesgue 
measure in M. d . Indeed, there exists an index i such that Lebesgue- almost every 
point x € R d belongs to L a v *^ (7a), i-e. satisfies 

dist (a;, < T Q ^(|a„| 1/Q ) for i.m. (k,n) e Z x M- 

Moreover, letting t/j denote the line spanned by the i-th vector of the canonical 
basis of M. d , we see that the hyperplanes H^ n are such that 

sup |{a; € f7, I dist(x, Hk , n ) < 1}| < 00. 

fcez 

We may therefore apply Theorem 3.6 in [35] and deduce that L o (o0 G G d_1,¥! (K ). 
To finish the proof, it suffices to consider a function 95 G $ Q with s v — a/^ a , and 
to recall that G d ~ 1 ' v (R d ) C c^- 1 ^**. Such functions exist; as a matter of fact, one 
may take <p(r) — T Q </j*(r) log(l/T Q ^(r)) where is given by LemmaO Q 

The above results being established, we are now in position to prove Theorems [2] 
and [3] This is the purpose of the last part of this section. 

11.3. End of the proof. Recall that we have to establish the following properties: 
The iso- Holder sets Ef(h) and the singularity sets E'^{h) of the Davenport series / 
have Hausdorff dimension equal to d—l + h/"/ a in every nonempty open set. We also 

need to show that the latter sets belong to the classes Q when h is positive. 

By virtue of Theorem D in [26] , this implies that the singularity sets have packing 
dimension equal to d in every nonempty open set, a feature that is also mentioned 
in the statement of Theorem [3] In view of various remarks written above, it only 
remains to consider the case where h < 7 and to establish the following three 
propositions. 

Proposition 9. For any real number h € [0,7 o ), 

max{dim H Et(h), dim H E' f (h)} < d — 1 + — . 

1 la 

Proof. We begin by making use of Lemma [4] The two inclusions ([43]) , combined 
with the fact that the sets Dj^ and Dm(o) are countable unions of hyperplanes, 
imply that 

dinin Ef(h) < max{d — 1, dimn E'Jh)} 
dimH E'f(h) < max < d — 1, inf dimu L a (a) 

1 a>h 
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To conclude, it suffices to apply Lemma [5] which, along with the alternate ex- 
pression ([20|) of 7 a , ensures that the dimension of L a (a) is bounded above by 
d-l + a/7 . □ 

In the next statement, G^ 1:V (R d ) denotes the extended class of sets with large 
intersection that is defined in terms of the initial class G d ~ 1 ' v (M d ) by the following 
condition: For all E C R d , 

E G G d_1 ' ¥> (R d ) 3E' G G d-1 '*'(R d ) £'c& 

The purpose of this extension is to avoid checking that the singularity sets are Gg- 
sets, which is inessential here. It is easy to see that the extended class G ' v (R d ) 
contains the initial class G d " 1 ' ip (R d ) and coincides with the latter on the Ga-sets, 
in view of [2lJ Proposition 1(e)]. Moreover, the extended class enjoys the same 
remarkable properties as the initial class: G 1:¥, (R d ) is closed under countable 
intersections and bi-Lipschitz mappings, and its members satisfy (|48j) . Moreover, 
when d — 1 + s v is positive, the class G 1 ' v3 (R d ) is included in the corresponding 
extended version Q d 1+s,p of Falconer's class. 

Proposition 10. Let us consider a real number h G (0,7 a ). Then, 

E' f (h)G p| G d - 1 - ip (R d )cg d - 1+hh \ 

'■pe<s> h 

Proof. Lemma [7] ensures that the set L a (h) belongs to the classes G d ~ 1 ' v (M. d ) asso- 
ciated with the functions (p G The same property holds for the set R d \ Dj/ a y, 
indeed, being the complement of a countable union of hyperplanes, this set is a 
G^-set with full Lebesgue measure in R d , and such a set belongs to all the classes 
G 9 (R d ), see [21] Proposition 11]. Using the stability under intersection of the classes 
of sets with large intersection, we deduce that 

La(h)\D J{a) e p| G^^jcs^ 1 *, 

where the last inclusion also follows from Lemma[7] To conclude, it suffices to make 
use of LemmaHJ which ensures that L a (h) \ -Dj(o) is a subset of E'^(h), see the first 
inclusion in (l44l). □ 



Our last statement gives a lower bound on the Hausdorff dimension of the sets 
Ef(h) and E'j{h) in every nonempty open subset of R d . Recall that Ef(0) contains 
the set -Dj( a ), by virtue of Lemma |U As 7 a is finite, the latter set is a dense 
countable union of hyperplanes, thereby having Hausdorff dimension at least d — 1 
in every nonempty open set. We may therefore restrict our attention to the positive 
values of h. 

Proposition 11. Let us consider a real number h G (0, 7 a ) and a nonempty open 
subset W ofR d . Then, 

mm{dim H (£/(/i) n W),dim H (E' f (h) n W)} >d-l + —. 

7a 

Proof. It follows from Proposition [TU] that the singularity set E'j(h) satisfies P5|) 
for any function ip G $ such that ip -< tp for some <p G Moreover, choosing a 
function ip for which s$ = h/j a , we also deduce from Lemma [5] that all the sets 
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L a (a), for a < h, have a vanishing H d '"^-rnass. Such a function i/j exists: It suffices 
to take i>{r) = T h ip*{r) (log(l/T ft ^(r))) 2 , and also y>(r) = T hV *(r) log(l/T h ^(r)), 
where is given by Lemma |6] We finally get 

n &-\rt {Ef(h) nw)> n A ~ x ^ (E' f (h) n w) = oo, 

where the first inequality is due to the second inclusion in (|44p . which appears in 
the statement of Lemma |U The result follows. □ 
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